Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 8 August 2011 (MN M£X style file v2.2) 



Self-similar scaling and evolution in the galaxy cluster 
X-ray Luminosity- Temperature relation. 

B. J. Maughan, 1 * P. A. Giles, 1 S. W. Randall, 2 C. Jones, 2 and W. R. Forman 2 

1 H. H. Wills Physics Laboratory, University of Bristol, Tyndall Ave, Bristol BS8 1TL, UK. 
2 Harvard-Smithsonian Center for Astrophysics, 60 Garden St, Cambridge, MA 02140, USA. 



o 

(N 
< 



O 

u 
6 



> 
o 
o 

(N 
00 

o 



X 



8 August 2011 



ABSTRACT 

Wc investigate the form and evolution of the X-ray luminosity-temperature (Lx — kT) 
relation of a sample of 114 galaxy clusters observed with Chandra at 0.1 < z < 1.3. 
The clusters were divided into subsamples based on their X-ray morphology or whether 
they host strong cool cores. We find that when the core regions are excluded, the most 
relaxed clusters (or those with the strongest cool cores) follow an Lx — kT relation with 
a slope that agrees well with simple self-similar expectations. This is supported by an 
analysis of the gas density profiles of the systems, which shows self-similar behaviour of 
the gas profiles of the relaxed clusters outside the core regions. By comparing our data 
with clusters in the REXCESS sample, which extends to lower masses, we find evidence 
that the self-similar behaviour of even the most relaxed clusters breaks at around 
3.5 keV. By contrast, the Lx — kT slopes of the subsamples of unrelaxed systems 
(or those without strong cool cores) are significantly steeper than the self-similar 
model, with lower mass systems appearing less luminous and higher mass systems 
appearing more luminous than the self-similar relation. We argue that these results 
are consistent with a model of non-gravitational energy input in clusters that combines 
central heating with entropy enhancements from merger shocks. Such enhancements 
could extend the impact of central energy input to larger radii in unrelaxed clusters, 
as suggested by our data. We also examine the evolution of the Lx — kT relation, and 
find that while the data appear inconsistent with simple self-similar evolution, the 
differences can be plausibly explained by selection bias, and thus we find no reason 
to rule out self-similar evolution. We show that the fraction of cool core clusters in 
our (non-representative) sample decreases at z > 0.5 and discuss the effect of this on 
measurements of the evolution in the Lx — kT relation. 

Key words: cosmology: observations - galaxies: clusters: general - galaxies: high- 
rcdshift - intcrgalactic medium - X-rays: galaxies 



1 INTRODUCTION 

The study of galaxy clusters is often motivated by one of 
two goals. First, clusters have been shown to be excellent 
cosmological probes, providing complementary and com- 
petitive constraints on cosmological pa rameters to those 
obtained using other techniques (e.g. I Allen et"aH 120081 : 
IVikhlinin et al.ll2009bl ). Second, galaxy clusters are unique 
laboratories in which to study extreme physical processes, 
such as cluster mergers or the interaction between rela- 
tivitic jets from active galactic nuclei (AGN) and the X- 
ray emitting plasma of the intra-cluster medium (ICM). 
In fact these goals are not exclusive. Accurate cluster 
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mass measurements are required for cosmological studies, 
but these must be estimated from observable properties 
which can be strongly affec ted by the non-gravitational pro- 
cesses at play in clusters dRandall. Sarazin fc Ricker|[ 2002: 



Rowley. Thomas fc Kavll2004lKravtsov. Vikhlinin fc Naeail 
20061 : iMagliocchetti fc Briiggenll2007l : lHartlev et alj|200a i. 



The scaling relations between different cluster observ- 
ables and between those observables and cluster masses are 
a case in point. The relationships between readily observable 
cluster properties such as X-ray luminosity (Lx) or temper- 
ature (kT) and mass can often be approximated as simple 
power laws, and can be used to provide "cheap" mass es- 
timates for larg e numbers of clusters out to high redshifts 
(|Maughanl 120071 ) . Power law scaling relations are expected 
under simplified models in which clusters are self-similar ob- 
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jects, having formed in single monolithic gravitational col- 
lapses and whose ICM is heated only by the shocks asso- 
ciated with the collapse ijKaiserl 1 198m . In this self-similar 
model, galaxy clusters and groups of all masses are identi- 
cal objects when scaled by their mas s. This is referred to 
as strong self-similarity (|Bower|[T997l ). and sets the power 
law slopes of the scaling relations which are not predicted 
to evolve with redshift. Evolution is expected, however, in 
the normalisation of the scaling relations, and is due (in the 
self-similar model) s olely to the changing den sity of the Uni- 
verse with redshift l)Brvan fc Normar] ll998). This redshift- 
dependent evolution in the normalisation of the scaling re- 
lations is referred to as weak self-similarity. 

The Lx — kT rel ation is the most well studied of these 
scaling relations Ce.g.lMitchell et al.lll979l;lMushotzkvl | l984 |; 
Edge fe Steward 1 199 it lMarkevitchlll998l ; lEttori et allhooi 



Pratt et al 



2009: 



Mittal et al.l 120111 ). as these two prop- 



erties can easily be measured directly and essentially in- 
dependently from X-ray data. Indeed, the observational 
consensus is that the slope of the Lx — kT relation is 
~ 3, significantly steeper than the self-similar prediction 
of 2, and appears to steepen further at galaxy group 
scale s (e.g. iHelsdon fc Ponmanl |2000| ; lOsmond fc Ponmanl 
|2004 ). The slope of the Lx — kT relation is the archety- 
pal example of the departure of clusters from the self- 
similar predictions. Other evidence for similarity break- 
ing is seen in the radial distribution of the ICM. Using 
simple models, iNeumann fc Arnaudl (|l999l ') found evidence 
for self-similarity of cluster surface brightness profiles out- 
side t he core regions, and lArnaud. Aghanim fc Neumann! 
J2002h extended this result to hi gher redshift clusters. 
IPonman. Cannon fc Navarro! l| 19991 ) showed that surface 
brightness profiles of the hotter (> 4 keV) clusters in their 
sample of relaxed systems were self similar outside of the 
core regions. Recently, with higher- quality data and more so- 
phisticated modellin g, the picture has become more detailed; 
ICroston et al.l (120081 ) used XMM-Newton observations to ex- 
amine the gas density profiles of a representative sample of 
31 local clusters and found that outside the core regions 
there is a temperature dependence of the gas density pro- 
files, but at larger radii, the profiles become self-similar for 
the most massive clu sters. Using the same sample of clusters, 
|Arnaud et alj|20ich also showed that the pressure profiles 
of the ICM have a universal, self-similar form, with low dis- 
persion outside the core regions. 

Perhaps the most detailed information on the processes 
responsible for similarity breaking has come from observa- 
tions of the entropy distribution of the ICM. Entropy is use- 
ful this regard because it is conserved in adiabatic processes, 
and so provides a good indicator of the thermal history of 
the gas. A key observation has been that the entropy pro- 
files of galaxy groups and clusters show departures from a 
self-similar power law, with profiles flattening in the core 
regions, indicating an a bsence of the expected low-entrop; 
gas in cluster cores (e.g. I Ponman. Cannon fc Navarrdll99' 



IVoit. Kav fc Brvanl 120051 ; IPratt et al.l 120101) . As this low- 
entropy gas is also the most luminous part of the ICM, a 
mechanism that removed or raised the entropy of this gas, 
and that had a proportionately larger effect in low mass 
systems, would also in principal explain the steeper than 
self-similar slope of the Lx — kT relation, and the similarity 
breaking in ICM density profiles. 



A great deal of work has been performed in investigat- 
ing the different processes that could be responsible for this 
simil arity breakin g, with cooling out of the low entropy gas 
(e.g. iBrvanl |2000| ) . preheating (e nergy input at early times 
in the formation of clusters; e.g. Bor gani et al.l 12002). and 
energy input from supernovae-driven galactic outflows all 
considered. A consensus is now emerging that energy input 
from AGN at high redshift, when they are in quasar-mode 
with high accretion rates, are responsible for providing a 
roughly constant (with halo mass) level of entropy injec- 
tion to the ICM, leadin g to the observed ICM structure 
and scaling relations (e . g.lBower. McCarthy fc Bensonll2008l ; 
IShort fc Thomasll2009l ; iMcCarthv et al.ll201ll ).~ 

Any study of the Lx — kT relation must include a con- 
sideration of the phenomenon of cool cores in clusters. These 
dense, cool regions in the centres of many clusters radiate 
efficiently in X-rays, and were long thought to be the sites 
of deposition of large quantities (h undreds of so lar masses 
per year) of condensed ICM (see iFabianl [l994, for a re- 
view). High resolution X-ray spectroscopy subsequently re- 
vealed that gas was not cooling out of the X-ray emitting 
phase in the large quantities predicted, with the bulk of 
it ins tead stabilising at around 1/3 of the global tempera- 
ture dPeterson et alj|200ll ; iKaastra et"afll200ll ; iDavid et all 
l200ll ). This effectively substituted the problem of explaining 
the fate of the cooling gas with the need for a mechanism 
to balance the cooling process to give the observed, lower 
cooling rates. This problem has been vigorously investigated 
over the past decade, with energy in put from AGN emerg- 
ing a s the most likely candidate (e.g. iMcNamara fc Nulsenl 
2007). This is a milder form of AGN input than the quasar- 
mode heating that drives similarity breaking. Here the en- 
ergy is thought to be input to the ICM via the bubbles 
inflated by the AGN jets and associated weak shocks, with 
the energy input coupled to the cooling rate by some feed- 
back mechanism to maintain a rough balance between heat- 
ing and cooling. Cool cores also have an important obser- 
vational effect in that they significantly increase (decrease) 
the measured global Lx (kT) compared to clusters without 
cool cores. The presence and strength of cool cores can thus 
influence the detectability of clusters in X-ray surveys, and 
complicate measurements of the form and evolution of the 
Lx — kT relation (and are a dominant source of scatter in 
luminosity scaling relations). 

Knowledge of the evolution of the Lx — kT relation is 
important as a step in the process of providing mass esti- 
mates for high-z clusters, understanding the selection func- 
tions of X-ray cluster surveys, and also for probing the his- 
tory of the heating mechanisms in clusters. For example, 
simulations have shown that AGN heating and preheating 
models give divergent predictions for the evolution of the 
Lx — kT relatio n, allowing the m odels to be distinguished 
observationally (|Short et al.ll201ol ) . Observational results are 
somewhat mixed, with some studies finding that evolution 
to z ~ 1 is consistent with self-similar predictions in which 
the evolution is driven by the increasing density of the 
Uni verse with redshift, leading to denser collapsed objects 



e.g.lVikhlinin et al-hOollLumb et al.l2004lMaughan et all 



2006), while other studies have found evidence for depar- 
tures from self-simila r evolution (e.g. lEttori et all 12004 ; 
iBranchesi et alj|2007l ). Recent work, however, has demon- 
strated the importance of considering selection bias in stud- 
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ies of the evolution of Lx scaling relations, showing that 
such biases can mimic or reduce sensitivity to departures 
from self-similar evolution Istanek et al.ll2006l : IPacaud et alJ 
120071 : iNord et alJliooj ; iMantz et alj|2010a[ ). 

In this paper we investigate the Lx — kT relation of 
a sample of 114 clusters covering wide temperature (2 < 
kT < 16 keV) and redshift (0.1 < z < 1.3) baselines to ex- 
amine the strong and weak self similarity of the cluster pop- 
ulatio n. The sample was first presented in iMaughan et alJ 
(2008, hereafter M08), wherein the analysis methods were 
described, and the evolution of the structural properties and 
metal abundance were investigated. Subsequently the sam- 
ple was used to show that Lx is a more precise mass est i- 
mator than had previously been thought (lMaughanll2007l ). 
The current paper is organised as follows: in §2 we re- 
view the sample and describ e updates to the analy sis and 
calibration of the data since IMaughan et al.l (|2008l ); in §3 
we discuss the classification of clusters based on dynami- 
cal and cool core state; in §4 and §5 respectively, strong 
and weak self-similarity are investigated, and in §6 and 
§7 the results are discussed and our conclusions are sum- 
marised. A ACDM cosmology of Ho = 70 km s -1 Mpc -1 (= 
WOh km s" 1 Mpc" 1 , and Q M = 0.3 (fi A = 0.7) is adopted 
throughout and all errors are quoted at the 68% confidence 
level. 



2 SAMPLE AND DATA REDUCTION 

The sample used for this study is a set of 114 galaxy clus- 
ters covering the redshift range 0.1 < z < 1.3 observed with 
Chandra ACIS-I, and originally presented in M08. In this 
work, the data have been reanalysed using updated versions 
of the CI AO software package (version 4.2), and the Chan- 
dra Calibration Database (version 4.3.0). The full analysis 
procedure is as described in M08, but the key steps are sum- 
marised in the following, along with any differences in pro- 
cedure from that work. 

The data were reprocessed from the level 1 events and 
cleaned and filtered. Blank-sky background files appropri- 
ate for each observation were also prepared and normalised 
to match the 9.5 — 12 keV count rate in the cluster data. 
Point sources were detected and excluded from all further 
analysis. The X-ray centroid of the cluster was determined, 
and radial profiles of the cluster and background data were 
used to determine the extent of the cluster emission, and to 
define local background regions free from cluster emission. 
Spectra were extracted from the cluster and blank-sky files 
in these local background regions and used to measure the 
difference in soft Galactic emission between the cluster and 
blank-s ky data. This differ ence spectrum was fit with a soft 
APEC IjSmith et al,ll2001f ) plasma model and included as 
a fixed component in the spectral fits to the cluster spec- 
tra, with the norma lisation scaled by extraction area (see 
IVikhlinin et al.ll2005l . for details). Fits to the cluster spectra 
were performed in the 0.6 — 9 keV band with an absorbed 
APEC m odel, with the absorbing co lumn fixed at the Galac- 
tic value (jDickev fc Lockmanlll99(jh . 

The gas density profile of each cluster was determined 
by converting the observed surface brightness profile (mea- 
sured in the 0.7 — 2 keV band) into a projected emissivity 
profile, which was then modelled by projecting a density 



model along the li ne of sight (see M08 for d etails). The model 
used was that of IVikhlinin et all (|2006bl . see that work for 
definitions of the parameters); 



nl(r/r c y 



x {l + r</r?)- h , 



(1) 



V ( 1 + r 2/ r 2)3/3- a /2 

a modification of the widely used /3- model, with added flexi- 
bility to fit a power law cusp in the core and a change in slope 
at large radii. Gas masses were then determined from Monte 
Carlo realisations of the projected emissivity profile based 
on the best fitting projected model to the original data. At 
each data point, a new randomised point was drawn from a 
Gaussian distribution centered on the model value at that 
point, with a standard deviation equal to the fractional mea- 
surement error on the original data point, multiplied by the 
model value. Note that this represents a minor change from 
M08 in which we randomised the original data to derive the 
uncertainties on the model parameters. The new approach is 
considered superior as we avoid adding noise to an already 
noisy observed profile. 

The cluster temperature, gas mass and -R500 (the radius 
enclosing a mean density of 500 times the critical density 
at the cluster's redshift) were then determined iteratively. 
The procedure followed was to extract a spectrum from 
within an estimated -R500 (with the central 15% of that ra- 
dius excluded) , integrate the gas density profile to determine 
the gas mass within the estimated R500, and thus calculate 
Fx (the product of k T and the gas mass, and a low sca tter 
proxy for total mass iKravtsov. Vikhlinin fc Nagaill2006l 1. A 
new value of -R500 was then estima ted from the Yx-M scaling 
relation of IVikhlinin et al.l i|2009ah . 



Ml, 



= E 



-2/5 



( Z )Ay 



Y x 



3 x 10 14 Mo keV 



(2) 



with A YM = 5.77 x 10 14 /i 1/2 M Q and B Y m = 0.57. 
Here, M500 is the mass within -R500 (allowing 
-R500 to be trivially computed), and E(z) = 
V^Ri + z) 3 + (1 - Q M - n A )(l + z) 2 + fi A , describ- 
ing the redshift evolution of the Hubble parameter. This 
Yx — M500 relation assumes self-similar evolution (as 
L -2 / 5 ), which has been sho wn to be a good description of 
observed clusters to z as 0.6 l|Maughanl 2007). Equation ([2]) 
is an updated version of the Yx-M relation used in M08, but 
the change is negligible here, leading to a < 2% increase in 
-R500 for the range of masses considered here. The process 
was repeated until -R500 converged. The temperature and 
luminosity were then measured from spectra extracted 
within -R500 both with and without the central 15% of that 
radius excluded. We use the notation Lx,tot and kT to t to 
indicate properties measured in the (0 — l)i?5oo aperture 
and Lx, c and kT c for those in the core-excised (0.15 — l)i?5oo 
aperture. 

The measured properties of the clusters are given in 
Table [T] All luminosities are bolometric. 



2.1 Comparison with M08 

In addition to the general methodological differences be- 
tween M08 and the current work that are described above, 
some improvements were made to the analyses of specific 
clusters. Specifically, for several clusters additional light- 
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Tabic 1: Summary of cluster properties. Lx tot an d kTtot were measured in 
the (0 — l)i?soo aperture and Lx,c and kT c were measured in the (0.15 — 
l)i?500 aperture. Clusters are in order of right ascension, and luminosities are 
bolometric. 



Cluster 


z 


R500 


fcTtot 


^X,tot 


kT c 








(Mpc) 


(kcV) 


(10 44 erg s" 1 ) 


(keV) 


(10 44 erg s- 1 ) 


MS0015. 9+1609 


0.541 


1 


.26 




50.6 ±0.6 


xnr — 

8-5-0.4 


35.5 ± 0.5 


RXJ0027.6+2616 


0.367 





.99 


_ „4- 1 'k 

O.Z_q y 


7.4 ±0.5 




5.6 ±0.5 


CLJ0030+2618 


0.500 





.84 


_i_n 


4.9 ±0.9 


-LI 7 

4-iii.o 


3.6 ±0.8 


A68 


0.255 


1 


.25 




17.7 ±0.5 


7.81 1 ;;; 


10.4 ±0.4 


A115 


0.197 


1 


.28 


°-°-o.i 


13.7±0.1 


D - '-0.3 


9.7±0.1 


A209 


0.206 


1 


.34 


_ in a 

'• z -0.4 


19.2 ±0.3 




13.1 ± 0.2 


CLJ0152.7-1357S 


0.831 


0. 


76 


A fi+ u -° 
4 -°-0.7 


9.1 ±0.6 


. 1 1 


7.1 ±0.6 


A267 


0.230 


1 


.07 


_i_n ^ 
4 ' y -0.3 


12.4 ±0.6 


7.3 ±0.5 


CLJ0152.7-1357N 


0.831 


0. 


,78 


_i_n 1 
c 1 +U. ( 

°- i -0.7 


11.8 ±0.6 


^■ y -0.8 


9.6 ±0.8 


MACSJ0159. 8-0849 


0.405 


1 


.32 


7 q + U - J 
'• y -0.3 


42.0 ±0.5 


10.2±g;f 


19.0 ±0.4 


CLJ0216-1747 


0.578 





78 


4-9 Q 

5.9IJ;? 


2.8 ±0.3 


I Q O 

°-°-i.s 


2.1 ±0.2 


RXJ0232.2-4420 


0.284 


1 


.33 


O.O_ 4 


30.7 ±0.9 


_|_ 1 A 

8-0±l;f 


16.7 ±0.7 


MACSJ0242. 5-2132 


0.314 


1 


.08 


, „_i_n 9 

4 ' D -0.2 


28.4 ±0.7 


_ __i_n 7 
-J-^_0.6 


7.8 ±0.3 


A383 


0.187 


1 


.03 


_i_n 1 
d - y -0.1 


9.9 ±0.2 


in ^ 
4 -°-0.3 


3.8 ±0.1 


MACSJ0257.6-2209 


0.322 


1 


17 


_ in ft 

'■ Z -0.5 


17.3 ±0.4 


ft 7+u.y 
°-'-0.6 


8.7 ±0.3 


MS0302. 7+1658 


0.424 





.81 


_i_n r 
d - 4 -0.4 


6.3 ±0.7 


q q-|-U.O 

°-°-0.6 


3.2 ±0.7 


CLJ0318-0302 


0.370 


0. 


.94 


-1-1 9 

°- 4 -0.9 


6.0 ±0.4 


4-1 7 
°- 4 -1.2 


3.9 ±0.3 


MACSJ0329.6-0211 


0.450 


1 


.01 


. --1-0 9 

4 -°-0.3 


28.6 ±0.8 


4 c+u.o 
4 -°-0.4 


11.6 ± 0.5 


MACSJ0404.6+1109 


0.355 


1 


09 


_i_n « 
°*°-0.5 


10.8 ±0.4 


_i_n 7 
°-°-0.5 


8.7 ±0.4 


MACSJ0429. 6-0253 


0.399 


1 


.08 


j_n 9 
°- z -0.2 


24.5 ±0.6 


in Q 

ft 

D.O_ y 


9.3 ±0.4 


RXJ0439.0+0715 


0.230 


1 


.14 


°" 4 -0.2 


16.1 ±0.3 


C- Q-|-U.4 

°-°-0.3 


8.7 ±0.3 


RXJ0439+0520 


0.208 





.97 


— +0.2 
^ ' -0.2 


9.3 ±0.4 


o.y_ 4 


3.3 ±0.3 


MACS J0451. 9+0006 


0.430 


1 


02 




16.9 ±0.8 


4-1 1 


11.1 ± 0.7 


A521 


0.253 


1 


.20 


-i_n 9 
4,y -0.2 


16.4 ±0.3 


4_n 9 

4.8tg;| 


13.8 ± 0.3 


A520 


0.199 


1 


.31 


_i_n 9 
D "°-0.2 


18.4 ±0.2 


4-fl ^ 

D -°-0.3 


14.5 ± 0.2 


MS0451. 6-0305 


0.550 


1 


.19 


6-8±g;l 


46.8 ±1.9 


4-1 9 

' -°-i.o 


31.3 ± 1.6 


CLJ0522-3625 


0.472 





82 


4 q + J-.u 
*-°-0.9 


3.3 ±0.3 


4- 1 A 

*-°-1.0 


2.5 ±0.2 


CLJ0542.8-4100 


0.642 





.91 


R 4+0.8 


11.3 ±0.5 


°' z -0.8 


8.6 ±0.4 


MACSJ0647.7+7015 


0.584 


1. 


26 


-I- 1 A 

io.9tJ:S 


43.8 ±1.1 


4-9 1 

AJ - d -1.6 


24.7 ±0.8 


1E0657-56 


0.296 


1 


.64 




79.3 ±0.6 


1 1 JA-V-O 

AJ -'-0.5 


55.4 ±0.5 


MACSJ0717.5+3745 


0.546 


1 


40 


1 1 o+ u - ' 
il ' Z -0.7 


79.3 ±1.2 


1 fi + 


59.0 ±0.9 


A586 


0.171 


1 


.23 


r> fj + 0.5 

'■ z -0.5 


15.2 ±0.4 


7 fi+ 0,8 
'■°-0.8 


7.6 ±0.4 


MACSJ0744.9+3927 


0.697 


1 


.08 


'■' -0.4 


48.9 ±0.9 


4_n 7 

8-i±8:6 


25.1 ± 0.7 


A665 


0.182 


1 


.42 


'■°-0.3 


23.4 ±0.3 


4_n zi 
7-8181 


16.6 ± 0.2 


A697 


0.282 


1 


49 


9-s±8:l 


40.4 ±0.6 


4-fl 8 

10.2lg-f 


26.4 ±0.5 


CLJ0848.7+4456 


0.574 





.57 


_ , _i_n 4 
^• 4 -0.3 


1.1 ± 0.2 


z - u -0.3 


0.9 ±0.2 


ZWCLJ1953 


0.320 


1 


.14 


_ „_i_n 
7 n ' ° 
'• u -0.6 


16.4 ± 0.4 


„ - 4-0 f\ 

°- A -0.6 


9.3 ±0.3 


CLJ0853+5759 


0.475 





.83 


-4-1 % 
°* U -0.9 


3.0 ±0.4 


4-1 ^ 

°- A -0.9 


2.6 ±0.4 


MS0906. 5+1110 


0.180 


1 


10 


j_n 9 
O.z_ 2 


9.0 ±0.2 


4-fl i 

4 -'-0.3 


5.0 ±0.1 


RXJ0910+5422 


1.110 





.48 


-1-1 zl 


2.8 ±0.3 


4- "1 Q 

z - ' -0.8 


2.2 ±0.4 


A773 


0.217 


1 


.30 


'■°-0.3 


18.3 ±0.3 


4-H zl 

7 4+ u - 4 
'■ 4 -0.4 


11.4 ±0.2 


A781 


0.298 


1 


.14 


c r + 0.6 

°-°-0.6 


11.8 ±0.5 


4-H 7 

°-°-0.5 


9.9 ±0.5 


CLJ0926+1242 


0.489 


0. 


84 


4 a+o.o 
4 - u -0.4 


4.6 ±0.3 


4-i n 

A r: + l-U 
4 "°-0.9 


3.2 ±0.3 


RBS797 


0.354 


1 


.19 


(\ 9+0.3 


50.3 ±0.9 


_ 11 n 
'' J -0.9 


13.0 ±0.6 


MACSJ0949. 8+1708 


0.384 


1 


.24 


in (S 

'■°-0.6 


28.6 ±0.7 


'• o -0.8 


17.4 ±0.6 


CLJ0956+4107 


0.587 





81 


1 0+O.8 
4 - 8 -0.7 


5.8 ± 0.3 


, n+0.7 
4 '°-0.5 


4.5 ± 0.4 


A907 


0.153 


1 


.14 


5 2+ ' 1 
°- z -0.1 


11.5±0.1 


°-*-0.2 


5.1±0.1 


MS1006.0+1202 


0.221 


1 


.12 


5 6+ ' 4 


7.7 ±0.2 


6-o±°:I 


5.1 ±0.2 


MS1008. 1-1224 


0.301 


1 


.04 


5 0+ ' 3 
°- u -0.3 


10.4 ±0.3 


4.6±°;t 


6.7 ±0.3 


ZW3146 


0.291 


1 


.28 


6 2+ 01 

°- z -o.i 


47.0 ±0.4 


7 «+ - 4 
'■8-0.4 


16.3 ± 0.3 


CLJ1113.1-2615 


0.725 


0. 


65 


q O+0.8 

d -°-0.6 


3.8 ±0.4 


3 6+ 13 


2.4 ±0.5 


A1204 


0.171 





.96 


34 +0.l 


9.9 ±0.2 


q 7 +0.3 
°-'-0.3 


2.7 ±0.1 


CLJ1117+1745 


0.548 


0. 


.71 


q 1+O.6 
^■-'■-O.B 


2.1 ±0.3 


q q+0.9 

^■■J-o.e 


1.8 ±0.2 


CLJ1120+4318 


0.600 





.93 


5 8 +0 ' 9 
°-°-0 .7 


13.4 ±0.8 


5 2+1-3 


8.7±0.7 


RXJ1121+2327 


0.562 


0. 


78 


q c+0.3 
— 0.3 


5.0 ±0.3 


q O+0.3 

°- z -0.3 


4.3 ±0.3 



continued on next page 



© 0000 RAS, MNRAS 000, 000-000 



The self-similar Lx — kT relation 5 



Table 1: continued 



Cluster 


z 


#500 


kT tot 




kT c 








(Mpc) 


(keV) 


(10 44 erg s- 1 ) 


(keV) 


(10 44 erg s- 1 ) 


A1240 


0.159 


0.92 


3-S-0.3 


1.8 + 0.1 


n c+U.3" 

d - 8 -0.3 


1.6 + 0.1 


MACSJ1131. 8-1955 


0.307 


1.43 


o 1+0.9 
8 - i -0.7 


31.5 + 1.0 


9 5+ 18 


21.9 + 0.9 


MS1137.5+6625 


0.782 


0.83 


6 6 +0 ' 8 
—0.7 


14.0 + 0.5 


6 5+ 1 - 2 
-°-1.0 


7.9 + 0.4 


MACSJ1149. 5+2223 


0.545 


1.28 


y - J -0.8 


47.2 + 1.2 


8.5_ ; 7 


35.3 ± 1.1 


A1413 


0.143 


1.30 


7 n+0.2 
'• z -0.2 


17.1 + 0.1 


7 i+0.3 

'■-■--0.3 


8.6 + 0.1 


CLJ1213+0253 


0.409 


0.79 


o r,+0.7 
°-°-0.6 


2.1 + 0.6 


o q+0.9 
J - y -0.8 


1.6 + 0.5 


RXJ1221+4918 


0.700 


0.90 


6 2 +0 ' 6 
D - z -0.6 


12.6 + 0.4 


5.9l°; 7 


10.1 ± 0.4 


CLJ1226.9+3332 


0.890 


0.97 


I0.3t£;i 


45.3 + 1.3 




24.5 ± 1.1 


RXJ1234.2+0947 


0.229 


1.15 


6.6± 2 ; 3 


6.2 + 0.4 


7 6+ 2 - 4 
'■ D -2.0 


5.5 + 0.3 


RDCS1252-29 


1.237 


0.55 


4 - 1-0.7 


6.5 + 0.7 


4.61°:? 


5.4 + 0.7 


A1682 


0.234 


1.20 


6 1+ 1 ' 3 


12.4 + 1.4 


5.8t 2 ;° 


9.7 + 1.2 


MACSJ1311. 0-0310 


0.494 


0.97 


= 5 +0.3 
O.O_ 2 


17.7 + 0.3 


6 5+ ' 8 
D -°-0.6 


7.2 + 0.2 


A1689 


0.183 


1.40 


o o+ - 3 
y - u -0.3 


39.4 + 0.3 


8 A+ 0A 
8 - 4 -0.3 


15.7 + 0.2 


RXJ1317.4+2911 


0.805 


0.51 




1.1 + 0.3 


9 9+O.8 
z - z -0.4 


0.9 + 0.7 


CLJ1334+5031 


0.620 


0.88 




7.6 + 0.9 




5.8 + 0.7 


A1763 


0.223 


1.38 


o 1+0.4 
°' 1 -0.4 


21.3 + 0.3 




14.5 ± 0.3 


RXJ1347.5-1145 


0.451 


1.51 


12 6 +0 - 4 
lz ' D -0.4 


155.6 ± 1.4 


14.2l| 4 


44.1 + 0.8 


RXJ1350.0+6007 


0.804 


0.72 


4 2+ ' 9 
*- z -0.7 


5.5 + 0.5 


4 -°-0.8 


4.6 + 0.5 


CLJ1354-0221 


0.546 


0.74 


, O+0.7 
d - 8 -0.6 


3.6 + 0.3 


o i+0.9 
6 - 1 -0.5 


3.0 + 0.3 


CLJ1415.1+3612 


1.030 


0.65 


= n+0.7 
°-°-0.6 


11.9 + 0.7 


4 o+O.O 

o q+0.5 
■3-y_0.4 


7.5 + 0.6 


RXJ1416+4446 


0.400 


0.87 


o pr+0.3 
°'°-0.2 


5.9 + 0.5 


3.6 + 0.4 


MACSJ1423. 8+2404 


0.543 


1.01 


6 2+°' 5 
D - z -0.4 


31.6 + 1.0 


6 +1 -° 


12.0 + 0.7 


A1914 


0.171 


1.40 


9 6+°' 3 
y - D -0.3 


34.3 + 0.3 


o cr+0.6 
°- b -0.4 


14.0 + 0.2 


A1942 


0.224 


0.97 


4 6+°' 3 


4.0 + 0.1 


4 9+0.3 
4 - z -0.3 


3.0 + 0.1 


MS1455.0+2232 


0.258 


1.07 


4 -°-o.i 


21.9 + 0.2 


4.7l°; 2 


6.4 + 0.1 


RXJ1504-0248 


0.215 


1.40 


7 1+0.2 


66.6 + 0.7 




15.1 ± 0.5 


A2034 


0.113 


1.23 


c c+0.1 

6 - 6 -0.2 


9.4 + 0.1 


6 3+ - 2 
O.O_ 2 


6.6 + 0.1 


A2069 


0.116 


1.21 


6 0+°' 2 

D.U_ 2 


6.4 + 0.1 


5 9+0.3 
°' y -0.3 


5.5 + 0.1 


RXJ1525+0958 


0.516 


0.84 


o 7 +0.3 
a -' -0.3 


6.7 + 0.4 


o c+0.3 
°-°-0.4 


5.7 + 0.4 


RXJ1532.9+3021 


0.345 


1.14 


5 1+0- 2 


38.0 + 0.9 


6 3+ 10 


11.7 + 0.5 


A2111 


0.229 


1.23 


6 4 +0 ' 6 
D - 4 -0.5 


11.4 + 0.3 


6 4 +0 - 7 


8.3 + 0.3 


A2125 


0.246 


0.79 


2 5+ 01 

z -°-o.i 


2.0 + 0.1 


2.4t°;» 


1.7 + 0.1 


A2163 


0.203 


1.91 


14 7+0-8 

i4 -'-0.9 


93.9 + 1.3 


15-2±1;1 


58.6 + 1.2 


MACSJ1621. 3+3810 


0.463 


1.03 


6 1+ ' 3 


19.8 + 0.6 


6 2 +0 - 5 
°- Z -0.5 


9.2 + 0.5 


MS1621.5+2640 


0.426 


1.04 


6 l+°' 6 


11.2 + 0.4 


°-°-0.6 


8.7 + 0.4 


A2204 


0.152 


1.44 


7 1+0.2 
'• i -0.2 


41.1 + 0.4 


8 4 +0 - 8 


13.1 ± 0.3 


A2218 


0.176 


1.21 


6 4 +0 ' 2 

-0.2 


13.3 + 0.1 


6 +0 - 3 


8.3 + 0.1 


CLJ1641+4001 


0.464 


0.77 


o c+0.4 
°-°-0.4 


2.7 + 0.3 


o c+0.6 

6 -°-o.5 


2.0 + 0.3 


RXJ1701+6414 


0.453 


0.88 


, O+0.3 
d ' 8 -0.3 


6.8 + 0.6 


4 1+0-5 


4.5 + 0.6 


RXJ1716.9+6708 


0.813 


0.81 


6 3+ 10 
D - o -0.8 


13.5 + 0.7 


7 +l.l 

°-'-0.9 


9.0 + 0.7 


A2259 


0.164 


1.11 


5 2+ ' 3 
°- z -0.3 


8.8 + 0.2 


r, 9+0.4 
°- z -0.4 


5.3 + 0.2 


RXJ1720.1+2638 


0.164 


1.27 


K q+01 

°- y -o.i 


22.3 + 0.2 


b ' 8 -0.3 


8.4 + 0.2 


MACSJ1720. 2+3536 


0.387 


1.15 


6 2+ ' 4 
D - z -0.3 


25.6 + 0.6 


7 9+O.9 
,,z -0.8 


11.5 ± 0.4 


A2261 


0.224 


1.34 


7 o+0.2 
' - J -0.2 


29.1 + 0.3 


7 o+0.4 
' -°-0.4 


14.3 ± 0.2 


A2294 


0.178 


1.35 


s Q+0.8 
o.»_ 7 


15.4 + 0.4 


8 4+1-1 

°- 4 -0.8 


9.4 + 0.3 


MACSJ1824.3+4309 


0.487 


0.90 


4 1+0-9 


5.7 + 0.7 


4 8+ 1 - 4 

4 -°-0.9 


4.5 + 0.5 


MACSJ1931. 8-2634 


0.352 


1.22 


5 9 +0 ' 3 

°- a -0.3 


46.8 + 0.9 


6 7+ 1 - 1 
"•'-0.7 


14.7 + 0.5 


RXJ2011. 3-5725 


0.279 


0.89 


4 +0 ' 2 
4 - u -0.2 


6.9 + 0.3 


3 6 +0 - 4 


3.0 + 0.2 


MS2053. 7-0449 


0.583 


0.76 


4 0+°' 4 


4.9 + 0.3 


3.91™ 


3.1 + 0.3 


MACSJ2129. 4-0741 


0.594 


1.12 


8 9+ 1 ' 1 


37.2 + 1.2 


S.3±\-{ 


22.3 ± 1.0 


RXJ2129.6+0005 


0.235 


1.22 


c: n + 0.2 
°- z -0.2 


21.4 + 0.5 


6 2 +0 - 6 
D - z -0.6 


9.8 + 0.4 


A2409 


0.148 


1.18 


c =+0.3 
°'°-0.2 


11.5 + 0.2 


= 7+0.4 
°- ' -0.4 


6.7 + 0.2 


MACSJ2228.5+2036 


0.412 


1.29 


8-2ta6 


33.4 + 0.8 




22.2 ± 0.7 


MACSJ2229. 7-2755 


0.324 


1.03 


4 1+0.2 


20.5 + 0.8 


°' u -0.7 


6.4 + 0.4 


MACSJ2245.0+2637 


0.301 


1.05 


r- 9+0.3 
°- z -0.3 


17.1 + 0.5 


4 9+0-5 

4 ' y -0.5 


6.8 + 0.4 


RXJ2247+0337 


0.200 


0.67 


9 7+0.7 
Z -'-0.5 


0.4 + 0.1 


o 9+0.9 

z - y -0.6 


0.3 + 0.1 


AS1063 


0.348 


1.56 


ii cr+0.6 
■ L1 -°-0.6 


93.6 + 1.1 


11 2+ 11 
1J - z -0.9 


42.1 ± 0.9 


CLJ2302.8+0844 


0.722 


0.76 


5 4+1.5 


5.2 + 0.4 


5 5+ 2 - 4 
0-0-1.4 


3.8 + 0.4 


A2631 


0.273 


1.26 


7 0+°' 7 

' - U -0.5 


19.5 + 0.5 


6 9+°' 8 
D - y -0.5 


14.1 ± 0.5 
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curve cleaning was performed by hand to remove periods of 
high background that were not detected in the M08 analysis. 
This generally led to a small decrease in cluster temperature, 
consistent with the removal of a harder spectral component. 
For two clusters (RXJ1701+6414, RXJ1525+0958) the de- 
crease was significant compared to the statistical uncertain- 
ties on kT (both decreased from « 5 keV to m 4 keV). For 
two other clusters the change in kT was large, but within 
the statistical uncertainties: CLJ0216-1747 decreased from 
fa 8 keV to « 6 keV and CLJ1334+5031 decreased from 
fa 16 keV to « 6 keV. For one cluster (CLJ1216+2633), the 
good time remaining after the new cleaning was too short 
for useful analysis, and this cluster was dropped from the 
sample, reducing the total number to 114 clusters. Finally, 
for AS1063, a redshift of 0.252 was erroneously used in M08, 
and this was corrected to 0.348 for the current analysis. 

There have been several significant updates to the 
Chandra calibration since the analysis presented in M08, 
which was based on CALDB version 3.2.3. The most signif- 
icant changes for the measurement of cluster temperatures 
and luminosities were updates to th e mirror effe c tive a rea, 
and the ACIS contamination model. iReese et all ((2010) re- 
cently examined the effect of Chandra calibration changes on 
the temperatures of galaxy clusters, finding that on average, 
temperatures measured with CALDB 3.1 were 6% higher 
than those measured with CALDB 4.2. The average temper- 
ature change (for kT measured within -R500) from M08 to the 
current analysis was < kT id/kT new >= 1.00 ± 0.06, where 
the uncertainty is the standard deviation of the clusters, and 
those clusters with significant individual changes discussed 
ab ove were exc l uded. This is smaller than the change found 
bv lReese etahl (|20ld l. but we note that unlike that work, 
our study was not optimised to study the effect of calibra- 
tion changes; there were numerous other changes between 
the M08 analysis and the current work, as discussed above. 

The bolometric luminosities measured within -R500 were 
also compared with M08, and the new luminosities were 
found to be higher by 8% on average (< L x ,oid/ Lx,new >= 
0.93±0.06). This is consistent with the « 9% decrease in the 
effective area at low energies introduced in CALDB 4.1.]Q. 

We note that all of the main results in the following sec- 
tions were also present with similar statistical significance in 
the data as presented in M08, so our results are not sensitive 
to the various calibration and analysis changes. 



3 QUANTIFYING DYNAMICAL STATE AND 
COOL CORE PRESENCE 

The dynamical state of the cluster was quantifie d by measur- 
ing th e centroid shift, following the method of I Poole et alj 
(2006). The centroid shift ((w)) was defined as the standard 
deviation of the distance between the X-ray peak and cen- 
troid, with the latter measured in a series of circular aper- 
tures centred on the X-ray peak and decreasing in steps of 
5% from -R500 to 0.05i?soo- For some clusters, -R500 extended 
beyond the edges of the detector in some directions. A max- 
imum radius (R w ) within which (w) could be measured was 



http : //asc .harvard. edu/ciao/why/caldb4 . 1 . l_hrma.html 
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Figure 1. Comparison of simple cool core proxies -F C ore (top 
panel) and cuspiness (bottom panel) with <r T . The dashed hori- 
zontal line in each plot indicates the value of each property used 
to divide the sample into CC and NCC clusters. The outlier at 
ctt ~ 8 is the merging, double cool core cluster A115. 

thus defined for each cluster, such that the maximum frac- 
tion of the area within this radius that was lost to excluded 
regions, or o ff the edge of the detector was 2%. Following 
iPoole et all |2006r i. the central 30 kpc were excluded from 
these centroid measurements to increase the sensitivity to 
faint structure. In an improvement to M08, the uncertain- 
ties on the centroid shifts were calculated from Monte Carlo 
randomisations of the X-ray images, with the measurement 
of (w) being repeated for 100 realisations of the input im- 
age with pixels randomised under a Poisson distribution. 
This takes into account the effect of the image noise, while 
in M08, the error was derived simply from the number of 
values used in the standard deviation calculation. 

In addition to quantifying cluster morphology with the 
centroid shift, two quantities were used to determine if a 
cluster was likely to host a cool core. In a recent comprehen- 
sive study, iHudson et all (|2010l ) tested 16 cool core probes, 
and concluded that for high quality data, a direct measure- 
ment of the central cooling time is preferred, but in lower 
quality data the cus piness of the gas density profile was rec- 
ommended (see also IVikhlinin et al.ll2006al ). As our sample 
spans a wide range of redshift and Lx, cuspiness is a more 
appropriate choice as it can be measured reliably for all clus- 
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Figure 2. Comparison of cluster dynamical state with cool core proxy. The centroid shift is plotted against <tt (left panel), -F core (centre 
panel) and cuspincss (right panel). The dashed line indicates the value of (w) = O.OO6R500 used to separate relaxed and unrelaxed 
clusters. 



ters. Cuspiness is defined as the logarithmic slope of the gas 
density profile at a radius of 0.04i?5oo, and was measured 
using the best fitting gas density models, with the uncer- 
tainty derived from the cuspiness measured for each of the 
Monte Carlo randomisations of the density profile. The sec- 
ond quantity used was the core flux ratio (Fcore), defined as 
the ratio of the bolometric unabsorbed flux from the cen- 
tral 0.15J?5oo to the total within i?soo- This quantity is the 
most robust to variations in the X-ray data quality and an- 
gular scale with redshift. F COIC values were also computed 
using rest-frame (0.5 — 2) keV band unabsorbed fluxes, but 
this made no significant difference, so bolometric fluxes were 
used. 

The most direct test of the presence of a cool core would 
be a temperature profile, but the data available are not suffi- 
cient for full temperature profiles for the majority of clusters 
in our sample. Instead, the projected temperature was mea- 
sured for each cluster in a core (0 — 0.15J?5oo) aperture and 
outer (0.15 — 0.3J?5oo) aperture, chosen to sample respec- 
tively the central dip and subsequent peak of the temper- 
ature profile of a typical cool core cluster (jVikhlinin et al.l 
2006b). The significance of the difference between the outer 
and core temperatures then gives a direct measurement of 
any cool core. This quantity (or) was defined as the differ- 
ence in temperatures divided by the quadrature sum of their 
measurement errors, such that positive ot indicates a cooler 
core region. This probe is sensitive to the data quality, and 
not just the strength of any cool core, so cool core clusters 
with low signal to noise observations could have ut ~ 0. 
However, large values of ctt are an unambiguous indicator 
of cool core presence and can be used to calibrate other cool 
core proxies which are less sensitive to data quality. The 
dynamical properties of the clusters are given in Table [2] 

Figure [T] shows the cool core proxies plotted against ot, 
and these plots were used to define the value of each proxy 
used to split the sample into clusters with and without cool 
cores (referred to as CC and NCC clusters hereafter) . Almost 
all clusters with ctt > 2 have F cola > 0.5 (i.e. more than half 
of the flux within -R500 comes from the central 15%), and 
we thus define F COIC — 0.5 as the border between CC and 
NCC clusters in the -Fcore parameter space, giving 31 CC and 
83 NCC clusters. Similarly, almost all clusters with ot > 2 
have cuspiness > 0.8, so this is defined as the CC/NCC 
boundary in cuspiness, giving 30 CC and 84 NCC clusters 



(this is simil ar to the value of 0.7 a dopted to define strong 
cool cores in IVikhlinin et~aT]|2006al ). 

There is substantial overlap between the subsamples 
defined in .Fcore and cuspiness, and unless specified other- 
wise, in the following sections all results are shown for the 
CC/NCC populations defined by F CO re. 

The dynamical state of the clusters is also related to 
their cool core status in the sense that cool core clusters 
tend to be dynamically relaxed. This is illustrated in figure 
[2] which shows the centroid shift parameter plotted against 
each of the cool core measures. The most relaxed clusters 
show evidence for cool cores in all of the different measures 
(albeit with a larger dispersion for cuspiness), and we adopt 
a threshold of (w) = O.OO6-R500 below which clusters are 
classed as "relaxed" (28 clusters), and above which as "un- 
relaxed" (86 clusters). This a minor update to th e value of 

w) = 0.005J?5oo used to segregate the clusters in iMaughanl 

2007h . 

We finally define a split into 21 clusters that passed the 
combined filter of all three of the above tests (F CO rc > 0.5, 
cuspiness > 0.8, (w) > O.OO6-R500), and 93 clusters that do 
not. These clusters that are both relaxed and pass both cool 
core filters are referred to as RCC clusters (for "relaxed cool 
core clusters"), while the complementary set are referred to 
as NRCC clusters. 



4 STRONG SELF-SIMILARITY 

The strong self-similarity of the Lx — kT relation of the 
population was investigated. Figure [3] shows the Lx — kT 
relation where both quantities were measured within -R500, 
including the core regions. A power law of the form 



L x = E(z)A LT (kT/T, 



(3) 



was fit to the data, with T* set at 6 keV (close to the 
median temperature of 5.9 keV for the full sample). The 
mode l was fit in log space usin g BCES orthogonal regres- 
sion l|Akritas fc Bershadvl Il996l ) , and the intrinsic scatter 
about the best fitting model was estimated by adding an 
extra error term in quadra ture to every da ta point until the 
reduced \ 2 was unity (see iMa ughan 20071. for details). The 
intrinsic scatter of a population is thus measured in the Lx 
direction about the best fitting Lx — kT relation for that 
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Tabic 2: Summary of the cluster dynamical properties, cr-p indicates the signif- 
icance of the temperature decrease in the core region, F core is the ratio of the 
bolomctric flux within 0.15i?5oo to that within -R50O1 cuspiness is the logarith- 
mic slope of the gas density profile at 0.04iJsoo, and (w) is the centroid shift, 
with R w giving the maximum radius within which (w) could be measured. 
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1 


,2 


0. 


212 


± 


0. 


022 


0. 


031 


■f 0. 

-0 


422 
024 


82.2 ± 10.2 


1 


00 


A267 


-1 


.3 





.415 


± 


0. 


036 





.105 


■f 0. 

-0 


041 
.040 


26.8 + 12.6 


1 


.00 


CLJ0152.7-1357N 


-0 


.4 





192 


± 


0. 


018 





.012 


-t-0 

-0 


029 
.009 


37.0 + 13.9 


1 


.00 


MACSJ0159.8-0849 


2 


.8 





.547 


± 


0. 


014 


1 


.147 


■f 0. 
-0 


011 

010 


2.8 + 0.3 


1 


.00 


CLJ0216-1747 








0. 


255 


± 


0. 


038 


0. 


073 


■f 1. 

-0 


109 
(KiS 


36.6 ± 12.7 


1 


00 


RXJ0232.2-4420 


1 


.1 





.455 


± 


0. 


027 





.984 


■f 0. 
-0 


024 
.038 


20.4 + 2.2 


1 


.00 


MACSJ0242.5-2132 


1 


1 


0. 


.726 


± 


0. 


065 


1 


.263 


+0 

-0 


047 

037 


6.0 + 1.6 


1 


.00 


A383 


4. 


.6 


0. 


.616 


± 


0. 


029 


1 


.221 


■f 0. 
-0 


()09 
.009 


1.9 + 0.3 


1 


.00 


MACSJ0257.6-2209 





.9 


0. 


.496 


± 


0. 


024 





531 


+0 

-0 


{171 
201 


7.6 + 1.9 


1 


00 


MS0302.7+1658 


-0 


.5 





.487 


± 


0. 


124 


0. 


.252 


-t-3. 
-0 


511 
.209 


17.3 ± 3.6 


1 


.00 


CLJ0318-0302 


-0 


.5 





339 


± 





034 





.062 


■f 3. 
-0 


450 
.049 


20.5 ± 7.7 


1 


.00 


MACSJ0329.6-0211 





.8 





.593 


± 


0. 


042 


1 


.213 


■f 0. 
-0 


(320 
.026 


14.0 + 2.6 


1 


00 


MACSJ0404.6+1109 


0. 


.6 





194 


± 





011 


0. 


.063 


+0 

-0 


(354 
056 


38.2 ± 7.5 


1 


00 


MACSJ0429.6-0253 


1 


.9 


0. 


.621 


± 


0. 


040 


1 


.252 


■f 0. 
-0 


016 


3.9 + 0.7 


1 


.00 


RXJ0439. 0+0715 





.4 





.460 


± 


0. 


020 


0. 


.769 


+0 

-0 


(317 
155 


8.0 + 2.4 


1 


.00 


RXJ0439+0520 


1 





0. 


.648 


± 


0. 


073 


1 


.240 


■f 0. 
-0 


(320 
.025 


5.7 + 1.6 


1 


.00 


MACSJ0451. 9+0006 


-1 


.3 


0. 


345 


± 


0. 


030 


0. 


.065 


+0 

-0 


.474 
053 


19.3 ± 8.0 


1 


00 


A521 





.1 





.154 


± 


0. 


004 





.559 


■f 0. 
-0 


003 
.094 


53.8 + 4.1 


1 


.00 


A520 


-0 


.6 


0. 


213 


± 





005 


0. 


.043 


+0 

-0 


(3(32 

.006 


63.6 ± 38.7 


1 


.00 


MS0451.6-0305 


-0 


.2 





.332 


± 


0. 


023 





233 


■f 0. 
-0 


137 

.203 


22.7 + 15.0 


1 


.00 


CLJ0522-3625 





3 


0. 


253 


± 


0. 


033 


0. 


071 


■f 0. 

-0 


702 
062 


12.8 ± 6.8 


1 


00 


CLJ0542.8-4100 


-0 


.9 





.242 


± 


0. 


016 


0. 


.059 


■f 0. 

-0 


178 
.051 


43.4 + 15.8 


1 


.00 


MACSJ0647.7+7015 


-0 


.3 


0. 


.437 


± 


0. 


021 





.368 


-t-0 

-0 


100 

150 


6.2 + 2.9 


1 


.00 


1E0657-56 


1 


.2 





.301 


± 


0. 


004 





.146 


-t-o. 

-0 


(314 
.022 


17.5 + 2.7 


1 


.00 


MACSJ0717.5+3745 





,7 





.256 


± 


0. 


006 


0. 


211 


+0 

-0 


086 
(352 


25.5 ± 12.6 


1 


00 


A586 


-0 


.2 





.502 


± 


0. 


030 


0. 


.391 


■f 0. 
-0 


.106 


3.3 + 1.7 


0. 


.75 


MACSJ0744.9+3927 





.7 





487 


± 





019 


0. 


455 


+0 

-0 


(158 
.036 


16.0 + 1.1 


1 


.00 


A665 





.8 





.288 


± 


0. 


005 





.123 


■f 0. 
-0 


(305 
.034 


32.8 ± 6.8 


0. 


.75 


A697 





3 





347 


± 


0. 


009 


0. 


151 


+0 

-0 


(308 
037 


6.0 + 4.5 


0. 


95 


CLJ0848. 7+4456 





.0 





.212 


± 


0. 


057 





.055 


■f 0. 

-0 


568 
047 


28.0 + 10.4 


1 


.00 


ZWCLJ1953 





.1 


0. 


433 


± 





020 


0. 


181 


+0 

-0 


015 
.067 


17.7 + 3.6 


1 


.00 


CLJ0853+5759 


0. 





0. 


.120 


± 


0. 


022 





.008 


■f 0. 

-0 


(310 
001 


51.6 ± 24.9 


1 


.00 


MS0906.5+1110 


-2 


.3 





443 


± 


0. 


015 


0. 


.309 


+0 

-0 


027 
.036 


64.5 ± 2.0 


0. 


90 


RXJ0910+5422 





.0 





.196 


± 


0. 


042 


0. 


.343 


■f 0. 
-0 


615 
333 


25.8 ± 12.8 


1 


.00 


A773 


-1 


.3 


0. 


.378 


± 





011 





.195 


-t-0 
-0 


(306 
008 


7.3 + 2.1 


1 


00 


A781 


-0 


.3 


0. 


.154 


± 


0. 


009 





.005 


-t-0. 
-0 


(304 
002 


61.2 ± 17.7 


1 


00 


CLJ0926+1242 





.8 





303 


± 


0. 


032 


1 


.034 


+0 

-0 


.168 


2.6 + 1.1 


1 


00 


RBS797 


1 


.5 





.743 


± 


0. 


054 


1 


.099 


■f 0. 
-0 


069 
.061 


2.3 ± 0.7 


1 


00 


MACSJ0949.8+1708 





.6 


0. 


.393 


± 


0. 


018 


0. 


.406 


+0 

-0 


051 
.310 


8.1 + 4.3 


1 


.00 


CLJ0956+4107 








0. 


.219 


± 


0. 


020 





.013 


-1-0. 

-0 


006 
006 


31.6 ± 14.5 


1 


.00 


A907 


3 


.2 


0. 


555 


± 


0. 


012 


0. 


.934 


+0 

-0 


012 
.016 


12.9 ± 0.9 


0. 


95 


MS1006.0+1202 


1 


4 





.340 


± 


0. 


013 





192 


+0 

-0 


192 
101 


21.1 + 4.5 


1 


00 


MS1008.1-1224 





.1 





358 


± 





021 


0. 


.196 


+0 

-0 


144 
.139 


50.6 ± 12.1 


1 


.00 


ZW3146 


5 


1 


0. 


.654 


± 


0. 


016 





.984 


+0 

-0 


014 

.013 


3.5 + 0.4 


1 


.00 


CLJ1113. 1-2615 





.8 


0. 


355 


± 





077 


0. 


.306 


+0 

-0 


570 
.276 


12.9 ± 5.1 


1 


00 


A1204 


2 


,7 





.727 


± 


0. 


051 


1 


.186 


-1-0. 

-0 


015 

.022 


6.3 + 1.4 


1 


.00 


CLJ1117+1745 





.0 





.134 


± 


0. 


028 


0. 


.205 


+0 

-0 


481 
.199 


17.8 ± 10.4 


1 


.00 


CLJ1120+4318 


-0 


.4 





.352 


± 


0. 


036 





.366 


-1-0. 

-0 


082 
354 


9.4 + 11.0 


1 


.00 


RXJ1121+2327 


-0 


.5 





146 


± 


0. 


013 


0. 


009 


-1-0. 

-0 


006 
004 


56.7 + 28.6 


1. 


.00 



continued on next page 
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Table 2: continued 



Cluster 


0"x 




F 


:ore 




cuspiness 


<w)(10- 3 iJ BC ,o) 


R-iu (-R500) 


A1240 











091 


± 





005 





003±5J 


001 

.001 


21.6 ± 6.2 


1 


00 


MACSJ1131. 8-1955 


i 








306 


± 





016 





183l c 


088 
.142 


29.2 ± 3.7 


1 


00 


MS1137.5+6625 





1 





435 


± 





030 





426l|: 


.115 
.128 


8.2 + 3.1 


1 


00 


MACSJ1149.5+2223 


-0 


9 





251 


± 





010 





076l c 


.024 
.064 


16.4 + 12.3 


1 


00 


A1413 


-0 


3 





498 


± 





006 





713l c 


.017 
.023 


5.2 + 1.5 





85 


CLJ1213+0253 











258 


± 





118 





094±J: 


80!) 
.082 


19.2 ± 10.7 


1 


00 


RXJ1221+4918 





6 





195 


± 





010 





040±J= 


030 
.028 


11.2 ± 12.3 


1 


00 


CLJ1226.9+3332 


1 


5 





459 


± 





028 





46015; 


076 
.132 


9.5 + 3.1 


1 


00 


RXJ1234.2+0947 


1 


1 





110 


± 





009 





0091J: 


036 
.005 


74.3 ± 18.0 


1 


00 


RDCS1252-29 


-0 


8 





168 


± 





029 





014l c 


359 
.009 


7.4 + 3.4 


1 


00 


A1682 





2 





222 


± 





039 





0261J: 


024 
.017 


11.2 ± 5.1 





80 


MACSJ1311. 0-0310 


2 


4 





594 


± 





029 





565l c 


090 
.071 


2.2 + 0.8 


1 


00 


A1689 





7 





600 


± 





012 





827lj 


010 
.018 


2.2 + 0.8 


1 


00 


RXJ1317.4+2911 











174 


± 





120 





0691J 


092 
.066 


30.4 + 11.4 


1 


00 


CLJ1334+5031 











226 


± 





040 





089±5; 


.754 
.070 


76.3 ± 37.2 


1 


00 


A1763 


-0 


7 





319 


± 





008 





203l c 


005 
.127 


6.5 + 5.4 





90 


RXJ1347.5-1145 


3 


1 





717 


± 





023 





897±j 


007 
.014 


5.9 + 0.4 


1 


00 


RXJ1350.0+6007 


1 


1 





161 


± 





024 





034l c 


.646 
.028 


34.9 ± 9.2 


1 


00 


CLJ1354-0221 











157 


± 





021 





025±c 


.206 
.018 


27.6 ± 8.2 


1 


00 


CLJ1415. 1+3612 


-0 


9 





369 


± 





038 





907l c 


042 
.053 


13.0 + 3.9 


1 


00 


RXJ1416+4446 





3 





394 


± 





061 





8631J: 


.081 
.159 


10.0 ± 2.3 


1 


00 


MACSJ1423.8+2404 





5 





621 


± 





051 


1 


246 + c 


.019 
.033 


3.1 ± 1.5 


1 


00 


A1914 


-1 


7 





592 


± 





015 





13lj;£ 


.013 
.016 


12.5 ± 1.1 





65 


A1942 


-1 


2 





256 


± 





012 





28915; 


.154 
.183 


8.0 + 1.4 


1 


00 


MS1455.0+2232 


4 


3 





709 


± 





020 


1 


048l|: 


.013 
.011 


3.5 + 0.2 


1 


00 


RXJ1504-0248 


2 


6 





774 


± 





037 


1 


2121J; 


.135 
.002 


1.4 + 0.8 


1 


00 


A2034 


-0 


3 





302 


± 





004 





3391J: 


.006 
.078 


12.2 ± 9.3 





75 


A2069 





3 





147 


± 





003 





038l c 


.120 
.026 


10.8 + 4.8 





75 


RXJ1525+0958 


-0 


3 





162 


± 





015 





0081J: 


009 
.004 


27.7 + 20.1 


1 


00 


RXJ1532. 9+3021 


1 


2 





691 


± 





054 





87ll c 


053 
.039 


2.8 + 1.5 


1 


00 


A2111 





6 





274 


± 





014 





067±c 


.019 
.038 


6.7 + 2.9 


1 


00 


A2125 


-1 


8 





125 


± 





009 





007±c 


.008 
.005 


22.3 ± 2.6 


1 


00 


A2163 





9 





376 


± 





010 





199l c 


.010 
.032 


25.0 + 7.4 





60 


MACSJ1621. 3+3810 


2 


5 





538 


± 





038 





995l c 


•>.$•> 
.043 


2.2 + 0.6 


1 


00 


MS1621. 5+2640 


-0 


1 





226 


± 





013 





05ll c 


009 

.045 


32.1 ± 7.2 


1 


00 


A2204 


3 


7 





681 


± 





023 


1 


342l c 


024 
.011 


1.3 + 0.4 





75 


A2218 


-1 


9 





373 


± 





007 





125l c 


005 
.009 


23.0 + 7.5 


1 


00 


CLJ1641+4001 





9 





261 


± 





051 





035l c 


068 
.022 


6.5 + 2.4 


1 


00 


RXJ1701+6414 





7 





328 


± 





052 





684±^ 


064 
.276 


10.1 ± 2.7 


1 


00 


RXJ1716. 9+6708 


-1 


3 





335 


± 





033 





340±c 


.267 
.319 


7.8 + 3.0 


1 


00 


A2259 


-0 








394 


± 





018 





165l{: 


.014 
.122 


8.9 + 6.8 





85 


RXJ1720. 1+2638 


5 


1 





625 


± 





018 


1 


0841J: 


014 
.035 


1.1 + 0.3 


1 


00 


MACSJ1720.2+3536 


1 


8 





552 


± 





029 





983lc 


028 
.044 


2.4 + 0.6 


1 


00 


A2261 


1 


4 





509 


± 





Oil 





682l c 


.044 
.023 


8.5 + 0.8 


1 


00 


A2294 


-0 


7 





388 


± 





019 





4091J: 


140 
.289 


6.1 ± 1.5 





65 


MACSJ1824.3+4309 











210 


± 





036 


1 


1661J: 


.246 
.252 


48.6 ± 9.2 


1 


00 


MACSJ1931. 8-2634 


2 








686 


± 





040 


1 


219l c 


.010 
.019 


3.5 + 0.9 


1 


00 


RXJ2011. 3-5725 


1 


3 





561 


± 





051 





492l c 


.062 
.107 


2.8 + 0.8 


1 


00 


MS2053.7-0449 





4 





354 


± 





035 





513l c 


.127 
.237 


10.2 ± 3.1 


1 


00 


MACSJ2129.4-0741 


-1 


4 





401 


± 





024 





184l c 


018 
.117 


15.1 ± 6.9 


1 


00 


RXJ2129. 6+0005 


2 








544 


± 





028 


1 


02ll c 


.059 
.048 


5.2 + 1.4 


1 


00 


A2409 


2 


3 





418 


± 





016 





322l c 


.086 
.233 


9.0 + 2.9 





65 


MACSJ2228.5+2036 





6 





335 


± 





014 





1821J: 


008 
.130 


26.2 ± 5.2 


1 


00 


MACSJ2229.7-2755 


1 


3 





686 


± 





082 


1 


319l c 


.014 
.021 


2.5 + 1.0 


1 


00 


MACSJ2245.0+2637 


1 


1 





603 


± 





050 





874l[: 


069 
.052 


3.4 + 0.8 


1 


00 


RXJ2247+0337 


-0 


6 





284 


± 





100 





020lc 


245 
.013 


13.3 + 4.2 


1 


00 


AS1063 





7 





549 


± 





016 





27ll c 


022 
.017 


7.5 + 1.5 


1 


00 


CLJ2302.8+0844 





2 





281 


± 





038 





34ll c 


230 
.316 


25.7 + 9.2 


1 


00 


A2631 


-0 


9 





276 


± 





012 





0661° 


088 
.049 


8.9 + 4.6 


1 


00 
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Figure 3. Lx — kT relation for the cluster sample, separated 
into relaxed/unrelaxed subsamples (top panel) or CC/NCC sub- 
samples (bottom panel) . Lx and kT were measured within -R500 1 
including the core regions. Luminosities are corrected for self- 
similar evolution as indicated on the ordinate axis. The lines show 
the best fitting power laws determined from a BCES orthogonal 
regression fit (see text for details). 



Figure 4. Lx — kT relation for the cluster sample, separated into 
relaxed/unrelaxed subsamples (top panel) or CC/NCC subsam- 
ples (bottom panel) . Lx and kT were measured within R500 1 with 
the central 0.15/Jsoo excluded. Luminosities are corrected for self- 
similar evolution as indicated on the ordinate axis. The lines show 
the best fitting power laws determined from a BCES orthogonal 
regression fit (see text for details). 



population, and is denoted as cr L | T . 

The Lx — kT relation was fit for the full cluster pop- 
ulation, and the CC/NCC and relaxed/unrelaxed subsets, 
and the best fitting parameters and measured intrinsic dis- 
persion are given in Table [3] In all cases, the normalisation 
measured for the CC or relaxed clusters was significantly 
higher than that of the NCC or unrelaxed clusters. This is 
in line with expectations, and the effect is clear in the dis- 
tribution of the clusters in the Lx — kT plane. The intrinsic 
scatter of each of the relaxed/CC subsamples is larger than 
the unrelaxed/NCC subsamples, and the highest scatter is 
that of the population as a whole. This is consistent with 
the relaxed/CC clusters and unrelaxed/NCC clusters form- 
ing two fairly distinct populations, offset in the Lx — kT 
plane. 

In order to remove the strong effect of the cool cores, 
the Lx — kT relation was also derived for core-excised clus- 
ter properties (i.e. both Lx and kT measured in the [0.15 — 
l]i?50o aperture). The resulting Lx — kT relations are plot- 
ted in Figure [4] again for the whole population, and for re- 
laxed/unrelaxed and CC/NCC subsets. The relations were 
fit as before and the best fitting parameters are given in 
Table [4] In these core-excised relations, the Lx — kT nor- 
malisations for all subsamples and the overall population 



are in good agreement, indicating that the core excision 
has removed the offset between the relaxed/CC and unre- 
laxed/NCC populations. 

In comm on with previous s tudies of luminosity scaling 
relati ons (e.g. lMarkevitchlll998l ; iMaughanlfeoOTl ; iPratt et all 
2009), the intrinsic scatter of the whole population was sig- 
nificantly reduced by the exclusion of the cores, from 67% to 
30% in this case. Indeed, the intrinsic scatter of each of the 
subsamples was reduced. The typical reduction in dispersion 
for the unrelaxed/NCC clusters was from ~ 40% to ~ 30%, 
while for the relaxed/CC clusters the effect was stronger, 
with a typical reduction from ~ 50% to ~ 15% (albeit with 
larger uncertainties on the intrinsic scatter). Also of note 
is that the 21 RCC clusters (corresponding to the most re- 
laxed clusters with strongest cool cores) have no measure- 
able intrinsic scatter when the cores are removed, indicating 
a very regular population. The core-excluded Lx — kT rela- 
tion for subsamples separated by the combined filter (RCC 
and NRCC) is shown in Fig. [5] 

The removal of the core regions also has the effect of re- 
ducing the slope of the Lx — kT relation for each subsample, 
although the effect is not always significant. The strongest 
effect is for the population fit as a whole, and is most likely 
due to the offset Lx — kT relations of the CC and NCC 
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category 


niter 




alt (iv erg s j 


£>LT 


TLT (%) 


all 


none 


114 


12.94 ± 1.00 


3.63 ± 0.27 


67.2 


± 4.8 


CC 


cuspiness > 0.8 


30 


26.30 ± 3.68 


3.33 ± 0.60 


54.6 


± 4.6 


NCC 


cuspiness < 0.8 


84 


9.75 ± 0.55 


3.18 ± 0.20 


36.8 


± 4.4 


CC 


F COTC > 0.5 


31 


23.28 ± 2.27 


2.53 ± 0.44 


50.6 


± 5.6 


NCC 


Fcorc< 0.5 


83 


10.02 ± 0.56 


3.19 ± 0.20 


38.1 


± 4.4 


relaxed 


(w}< 0.006 


28 


22.23 ± 2.98 


3.10 ± 0.61 


61.1 


± 6.1 


urelaxed 


(w}> 0.006 


86 


10.54 ± 0.70 


3.26 ± 0.22 


48.5 


± 6.9 


RCC 


combined 


21 


27.35 ± 2.64 


2.44 ± 0.43 


38.0 


± 3.2 


NRCC 


combined 


93 


10.57 ± 0.64 


3.22 ± 0.20 


46.9 


± 6.0 



Table 3. Lx — kT relations for properties measured with core regions included ([0 — l]-Rsoo aperture). The first column gives the category 
of each subsample, and the second column gives the filter used to define that subsample. The combined filter refers to clusters with 
cuspiness> 0.8, -F C ore > 0.5, and (w) < 0.006. The third column gives the number of clusters in each subsample. The fourth, fifth and 
sixth columns give the normalisation and slope of the best fitting model and the intrinsic scatter of the data about that model. The 
model used was L x = E(z)A LT (kT/T,) B ^T % with T» = 6 kcV in all cases. 



category filter N A LT (10 44 erg s" 1 ) B LT cr LT (%) 



all 


none 


114 


6.98 ± 0.30 


2.72 ± 0.18 


29.7 ± 4.2 


CC 


cuspiness > 0.8 


30 


7.16 ± 0.33 


2.15 ± 0.17 


13.0 ± 11.7 


NCC 


cuspiness < 0.8 


84 


6.98 ± 0.39 


2.82 ± 0.21 


33.1 ± 5.3 


CC 


F corc > 0.5 


31 


7.05 ± 0.34 


2.10 ± 0.13 


19.2 ± 5.1 


NCC 


F corc < 0.5 


83 


7.18 ± 0.38 


2.86 ± 0.22 


31.4 ± 6.2 


relaxed 


(w)< 0.006 


28 


6.71 ± 0.34 


2.12 ± 0.17 


12.5 ± 6.5 


urelaxed 


(w)> 0.006 


86 


7.28 ± 0.39 


2.86 ± 0.21 


32.4 ± 5.9 


RCC 


combined 


21 


7.28 ± 0.26 


1.90 ± 0.14 


1.8 ± 3.1 


NRCC 


combined 


93 


7.10 ± 0.36 


2.82 ± 0.19 


32.5 ± 5.1 



Table 4. Lx — kT relations for properties measured with core regions excluded ([0.15 — l]i?soo aperture). The first column gives the 
category of each subsample, and the second column gives the filter used to define that subsample. The combined filter refers to clusters 
with cuspiness > 0.8, -F C orc > 0.5, and (w) < 0.006. The third column gives the number of clusters in each subsample. The fourth, fifth 
and sixth columns give the normalisation and slope of the best fitting model and the intrinsic scatter of the data about that model. The 
model used was L x = E(z)A LT (kT /T t ) BLT , with T* = 6 keV in all cases. 

clusters, combined with different distributions with temper- 
ature of the CC and NCC clusters. Specifically, the CC clus- 
ters which are offset above the rest of the population in the 
Lx — kT plane are generally hotter in our sample (the me- 
dian core-excised kT is 6.5 keV for CCs and 5.5 keV for 
NCCs). When the cores are excised, these clusters move sig- 
nificantly towards lower Lx and higher kT, giving rise to 
the flatter slope than when cores are not excised. 

Perhaps of greater interest is that with cores excised, 
a difference in Lx — kT slope is apparent between the re- 
laxed/CC and unrelaxed/NCC populations. In all cases, the 
best fitting slope to the relaxed/CC clusters is shallower 
than that of the complementary population (see Fig. [5]). 
The significance of the difference ranges from « 2.5a to 
~ 4a, with the strongest difference found for the populations 
split by the combined filter. It is striking that all of the re- 
laxed/CC subsamples have Lx — kT slopes consistent with 
the self-similar prediction of 2. This is a first indication that 
this part of the cluster population obeys self-similar scaling 
laws. This was investigated further via an examination of 
the cluster gas density profiles. 



4.1 Similarity of gas density profiles 

In order to investigate the apparent self-similarity of the 
Lx — kT relation for the relaxed/CC clusters, the structure 



of the gas density profiles were examined. The best fitting 
gas density models for each cluster were scaled by divid- 
ing the radial coordinate by -R500, and dividing the densities 
by E 2 (z) to remove the expected self-similar evolution (see 
e.g. ICroston et al.ll2008l . for examples of the effects of these 
scalings). After these scalings, self-similar profiles would be 
identical. The scaled profiles are plotted in Figure [U with 
the relaxed and unrelaxed clusters indicated. In both cases, 
the dispersion of the profiles is largest in the central regions, 
reducing at larger radii. For both subsamples, the dispersion 
in the profiles was significantly larger when the E 2 (z) scal- 
ing was not applied. It is apparent that the two subgroups 
have different gas density structures, with the relaxed clus- 
ters more centrally concentrated, and the unrelaxed clusters 
more diffuse. 

These results demonstrate a lack of similarity between 
the gas density profiles of clusters selected as relaxed and 
unrelaxed. Next, the degree of self similarity within those 
subgroups was investigated. Figure [7] shows the scaled gas 
density profiles of the relaxed and unrelaxed clusters, colour- 
coded by system temperature (measured with the core ex- 
cised). For the unrelaxed clusters, there is evidence for a 
temperature dependence of the profiles, with cooler systems 
tending to have lower densities than hotter systems at a 
given radius. For the relaxed clusters, no such trend is ap- 
parent. The temperature dependence of the density profiles 
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Figure 5. Lx — kT relation for the cluster sample, split using our combined filter into a subsample containing the most relaxed clusters 
with strongest cool cores (RCC), and the complementary subsample (NRCC). Lx an d kT were measured within -R50O1 with the central 
0.15-Rsoo excluded. Luminosities are corrected for self-similar evolution as indicated on the ordinate axis. The lines show the best fitting 
power laws determined from a BCES orthogonal regression fit (see text for details). 



was examined in more detail by measuring the scaled den- 
sity of each cluster at a fixed fraction of -R500, and plotting 
them against system temperature. An example plot for den- 
sities measured at 0.3-Rsoo is shown in Figure [8] Despite the 
different temperature distributions of the two populations, 
the unrelaxed clusters clearly have a stronger dependence of 
gas density on temperature than the relaxed clusters. 

The data for each subsample were then fit with a power 
law of the form: 

p(r/R 5O0 )E-' 2 (z) = (kT)\ (4) 

where we refer to r\ as the "similarity index" of the den- 
sity profiles (rj = corresponds to self-similar profiles). This 
procedure was repeated at different values of r to produce 
the similarity index profile shown in Figure(9] The similarity 
index profile shows that the gas density profiles of the unre- 
laxed clusters have a strong temperature dependence in the 
core regions that becomes weaker with radius, approaching 
self-similarity by ~ 0.7/Zsoo- A similar trend of decreasing 
simila rity index with radius was observed bv lCroston et all 
(2008) who found values of 0.5 and 0.25 at 0.37?soo and 
0.7-Rsoo respectively for their representative sample of clus- 
ters. 

In contrast, the similarity index profile of the relaxed 
clusters is quite flat, with no strong dependence on temper- 
ature at any radius. As a test of the sensitivity of this dif- 




Figure 6. Scaled gas density profiles of the cluster sample. The 
profiles are the best fitting models to the observed cluster emis- 
sivity profiles (see text for details), and are scaled in radius by 
-R500 and in density by E 2 (z). The solid lines indicate the unre- 
laxed clusters while the dashed lines are the profiles of the relaxed 
clusters. 



ference in similarity index to the absence of relaxed clusters 
with low temperatures, unrelaxed clusters with kT< 3.5 keV 
were excluded from the fits. This made no significant differ- 
ence to the slopes. 
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Figure 7. Scaled gas density profiles of the clusters, colour coded 
by system temperature (as measured in the [0.15 — l]i?500 aper- 
ture). The top panel plots the relaxed clusters and the bottom 
panel plots the unrelaxed clusters. 
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Figure 8. Scaled gas density at 0.3i?5oo is plotted against system 
temperature (as measured in the [0.15— l]i?500 aperture). Filled 
circles and hollow triangles indicate the unrelaxed and relaxed 
clusters respectively. 



Very similar results were found for density profiles when 
the CC and NCC samples were compared, although we con- 
sider the separation on (w) to be preferable as this is less 
directly linked to the radial structure of the gas density than 
either the Fcorc or cuspiness measurements. 



Figure 9. Similarity index (the slope of the power law depen- 
dence of scaled gas density on system temperature) is plotted 
as a function of scaled radius. Filled circles and hollow triangles 
indicate the unrelaxed and relaxed clusters respectively. The in- 
nermost points are slightly offset in radius for clarity. 



5 EVOLUTION OF THE L> 



kT RELATION 



Thus far, it has been assumed that the redshift evolution 
of the Lx — kT relation is self-similar, i.e. the slope is inde- 
pendent of redshift, while the normalisation increases with 
increasing redshift due to the increasing density of the clus- 
ters. As discussed in the introduction, there is some debate 
in the literature as to whether this evolution model is a good 
description of observed clusters, with evidence that selection 
biases play a significant role in evolution studies. 

To investigate this further, the evolution in the Lx — kT 
relation of the current Chandra cluster sample was exam- 
ined. In order to measure evolution, a low redshift baseline 
must be defined. While many exist in the literature, the 
limited fields of view of Chandra and XMM-Newton mean 
that at z < 0.1, the available Lx — kT relations were de- 
rived from earlier missions such as ROSAT and ASCA. For 
internal consistency, the Lx — kT relation of the clusters in 
the sample with 0.1 < z < 0.2 was fit as described above 
and used as the local baseline for our evolution study. In 
this fit, and all of the following evolution study, the core- 
excised properties were used. The parameters for the local 
L x - kT relation were Alt = (6.05 ± 0.04) x 10 44 erg s' 1 
and B LT = 2.51 ± 0.29. 

For each cluster, the ratio of its observed luminosity to 
that predicted by self similar evolution of the local relation 
was computed, and termed Al- Figure ITOl shows the plot 
of Al against redshift. In this plot, clusters following self- 
similar evolution would have Al = 1. Clusters more (less) 
luminous than predicted by the self-similar evolution of our 
local relation have Al > 1 (Al < 1). Figure [101 shows that 
relative to the local subset of clusters, those at intermediate 
redshifts (0.2 < z < 0.6) are significantly more luminous 
on average than the self-similar evolution model predicts. 
At high redshifts (z > 0.6) the clusters agree with the self- 
similar evolution model. 

Similar results were reported bv lBranchesi et all (|2007l . 
based in part on data from M08), and was suggested as 
evi dence against the self- similar model. However, as noted 
by iBranchesi et all (|2007l ). selection bias could be respon- 
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Figure 10. The ratio of the observed, core excised luminosity to 
that predicted by our local Lx — kT relation, including self similar 
evolution, is plotted as a function of redshift. Clusters evolving 
self-similar ly should scatter around Aj, = 1. Also shown are the 
weighted mean in several redshift bins, and the possible effect 
on Aj, due to Eddington bias in our simplistic bias simulations. 



sible for distort ing the observed evolution (see also §5.2 in 
lMaughan|[2007l ). 

At this point, we note that there is some ambiguity in 
the terms used to describe different selection biases in the 
literature on cluster scaling relations, and so we pause for a 
brief discussion of selection bias and terminology. 



5.1 Selection biases and terminology 

We consider two types of bias that affect flux-selected sam- 
ples of objects. Either or both of these biases are often 
referred to simply as "Malmquist bias" in the literature, 
while occa s ionall y the term "Eddington bias" is also used. 
iTeerikorpil (|l997l i gives a detailed overview of these selec- 
tion biases and advocates the use of the terms Malmquist 
bias of the first and second kind. For clarity, we use the 
term "Malmquist bias" to refer to the effect that arises due 
to more luminous objects being detectable to greater dis- 
tances, and "Eddington bias" to refer to the effect due to 
the differential scattering of objects across the threshold of a 
flux-limited sample. We will now discuss each of these biases 
in more detail in the following. 

Let us consider a population of objects whose luminosi- 
ties L axe correlated to their masses M, with some intrin- 
sic scatter in luminosity at a given mass (parametrised by 
0l\m)- A sample of these objects whose measured flux is 
greater than some flux limit Fu m is then defined. In the 
idealised case where there is no scatter, the Malmquist bias 
will lead to the most massive, luminous objects being over- 
represented in the sample, due to the larger search volume 
in which they may be detected. This is easily corrected for 
with knowledge of the survey selection function [the survey 
volume as a function of luminosity; V(L)] to weight the ob- 
jects by their search volume. In the more realistic case when 
there is scatter in L for fixed M, those objects with higher 
than average L for their mass, will be more numerous in the 
sample than those with lower than average L, again due to 
the larger survey volume. The mean luminosity of objects in 



the sample is thus biased higher than the population average 
and the size of the bias depends on oum- 

Turning now to the Eddington bias, let us consider ob- 
jects close to the flux limit (by "close" we mean objects 
whose scatter in L has a non-negligible chance of moving 
them across the flux limit). Near the flux limit there ex- 
ist objects whose mass corresponds to an average L that 
would be below the flux limit at their redshift, but whose 
observed L is higher than average, and so are included in 
the sample (we refer to these as overluminous). Similarly 
there exist objects whose mass corresponds to an average 
L that would be above the flux limit at their redshift, but 
whose observed L is lower than average, and so are not in- 
cluded in the sample (we refer to these as underluminous). 
This means that objects included in the sample are biased 
to above average L. So far, this is equivalent to a special 
case of the Malmquist bias, with the survey volume drop- 
ping to zero for objects whose observed luminosity is below 
the flux limit. What distinguishes the Eddington bias is an 
additional factor that arises when the number density of ob- 
jects is a decreasing function of M and L (as is the case 
for galaxy cluster mass and luminosity functions). In this 
case, there are larger numbers of lower mass, overluminous 
objects which are included in the sample than higher mass, 
underluminous objects which drop out of the sample. This 
gives a further bias towards above-average luminosity in the 
final sample. The size of the Eddington bias depends on 
CTh\Mi an d on the slope of the mass function at the mass 
whose average luminosity corresponds to the flux limit at 
the redshift in question. 

In the above discussion, we assumed that o L \ M was due 
solely to an intrinsic dispersion in L at fixed M. Similar 
biases occur if L is subject to measurement errors. In this 
case, even in the absence of intrinsic variation of L at fixed 
M, there is an observed variation in L which has the same 
effect as the biases described above. A complicating factor 
is that the size of o L \ M will typically depend on L due to 
the larger measurement errors on fainter systems. This in- 
creases the difference between the number of overluminous 
objects moving into the sample from below the flux limit and 
underluminous objects moving out of the sample, further in- 
creasing the bias. If the objects are reobserved the average 
L will be unbiased for fixed M, as the original L values were 
randomly drawn from the possible L distribution (this is in 
contrast to the case where the variation in L at fixed M is 
intrinsic). The number of objects in the sample will remain 
biased, however, unless objects below the flux limit are also 
reobserved, and the sample redefined. 



5.2 The effect of selection bias on cluster scaling 
relations 

In studies of the evolving mass function of clusters based 
on X-ray flux limited samples, both the Malmquist bias 
and Eddington bias must be accounted for. Typically, the 
Malmquist bias as defined here is included in modelling the 
survey selection function, while the Eddington bias is taken 
into account when modelling the mass function (see e.g . 
iMarkevitcbl 19981 ; IVikhlinin et al.ll2009al ; iMantz et al1l2010bl . 
for useful discussions of the effects of these biases). When 
X-ray selected clusters of galaxies are used for studying the 
scaling relation between Lx and another property (such as 
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kT), both the Malmquist and Eddington biases will enhance 
the mean Lx for a given kT (or mass). The effect of this can 
be to bias the normalisation of the Lx — kT relation high, 
and bias the slope low. The effect on the slope arises because 
the lower temperature systems are closer to the flux limit at 
a given redshift, and so will suffer stronger bias. This effect 
is lessened if a wide redshift range is considered, as higher 
temperature objects will also tend to be close to the flux 
limit. As described above, the size of the total bias in the 
Lx — kT relation will depend on (Tl\m> an d on the slope of 
the mass function at the flux limit of the survey. Note that 
while a L i M can be significantly reduced by excluding clus- 
ter cores for analysis purposes, the core-included Lx must 
be used in consideration of bias, as generally clusters are 
selected (or not) on the basis of their total Lx- 

In the sample presented here, the clusters are drawn 
from a heterogeneous archive, and do not represent a statis- 
tically complete sample. However, most of the clusters were 
originally detected in X-ray surveys, and are thus subject to 
Malmquist and Eddington bias. A rigorous treatment of the 
bias is impossible without a well-defined selection function, 
but some simple simulations were performed to investigate 
the plausible effect that selection bias could have on the ob- 
served evolution of the Lx — kT relation in this sample. The 
process used was to simulate a self-similar population of clus- 
ters, and apply selection functi ons representativ e of different 
cluster surveys. As discussed in lMaughanl (|2007l ). we approx- 
imate the selection of clusters in our sample by three very 
broad classes of X-ray surveys: at low redshift (z « 0.15) 
many clusters were selected in shallow all-sky X-ray sur- 
veys, for which we use the Brightest Cluster Survey (BCS; 
lEbeling et aTlll99St ) flux limit as a template (a bolomet- 
ric flux of 8.3 x 10 -12 erg s _1 cm -2 ); at moderate redshift 
(z ~ 0.45) most of the clust ers come from deeper surveys 
for which we use the MACS (|Ebeling. Edge fc Henrvll200ll . 
; bolometric flux limit of 3.4x 10 -12 erg s -1 cm -2 ) as a tem- 
plate (our sample contains 24 MACS clusters with a median 
redshift of 0.42); at high redshifts (z ~ 0.85) most clus- 
ters in our sample were detected in deep surveys of ROSAT 
pointings, for which we use the WARPS (|Scharf et al.ll 19971. 
bolometric flux limit of 9.1 x 10~ 14 erg s _1 cm -2 ) as a tem- 
plate. 

At each of the approximate redshifts corresponding to 
the above survey c lasses, a mass functi on was constructed 
using the model of I Jenkins et al.l l|200ll . see their Equation 
9). This model is based on fitting results from cosmological 
simulations, and is simil ar in form to the stand ard Press- 
Schechter mass function (Press & Schcchtcr 1974), although 
it provides a significantly improved fit to the simulations. 
The redshift range for each survey class was broken into 
smaller redshift bins, and clusters were drawn with random 
redshifts and masses (from the above mass function) until we 
obtained a set number of detected clusters within each bin 
(where a cluster is "detected" if its scattered total flux falls 
above the flux limit for the appropriate survey). The cluster 
luminosities were determined as follows. Each cluster was 
assigned a total (i.e. core included) luminosity Lx,tot and a 
core excis e d lum inosity Lx, c , using the Lx — M relations of 
iMaughanl l|2007h . Importantly, these relations include self 
similar evolution, so the luminosity at fixed temperature 
evolves smoothly as E(z). Lx,tot was then randomised un- 
der a lognormal distribution with a = 39% to give Lx,tot> 
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Figure 11. Plot showing the correlation between the ratio of 
observed to predicted Lx, when Lx and kT are measured with 
(A fj X tot ) and without ( A l x c ) the core regions included. 

reproducing the intrinsic dispersion of the IMaughanl (|2007l ) 
Lx,tot — M relation. The value of Lx, c was also shifted in 
the same direction by the same fractional amount but scaled 
by 17 /39 (the ratio of the intrinsic scatters in the IMaughanl 
120071 , Lx, c -M and Lx,tot-M relations) to give L' x c . This gives 
consistent values of the total and core excised luminosities 
for each cluster, including intrinsic scatter. By scaling Lx,c 
in this way, we are making the assumption that the dis- 
persion in Lx,tot is correlated with that in Lx, c - This is 
supported by Fig. 1111 which shows Az, xtot plotted against 
Al x c for the observed clusters (where A_t x tot is the ratio 
of the observed Lx,tot to that predicted by the low-redshift 
Lx,tot— kT relation, and A_l Xc is the same but for core- 
excised quantities), and shows a positive correlation. The 
weak correlation for the CC clusters in the plot is testa- 
ment to the efficacy of removing the core of those systems 
in reducing the dispersion in the Lx — kT relation. Finally, 
for each cluster, a tempera ture was assigned, derived from 
the IVikhlinm" et all (|2006bT ) M - kT relation, again assum- 
ing self similar evolution, but neglecting any intrinsic scat- 
ter. Thus, at each of our redshifts of interest, we have pro- 
duced a population of clusters whose observed properties 
(L x , t ot! ix,c an d kT) are (by construction) described per- 
fectly by a consistent set of self-similarly evolving scaling 
relations, with intrinsic dispersion in Lx. 

The flux limit of each of our three representative surveys 
was converted to a luminosity limit at the appropriate red- 
shift, and clusters with L x . to t below that limit were rejected. 
In this way, the effects of the selection biases are reproduced. 
For the clusters surviving the survey selections, the ratio of 
L XjC to Lx, c was calculated, which corresponds to the Al of 
the observed clusters. Figure [TD] shows the mean Al for the 
simulated clusters at each redshift. The departures from self 
similar evolution (Al = 1) are solely due to the selection 
biases, and as the intrinsic scatter o l \m was constant with 
redshift in our simulations, the variation in Al with redshift 
is simply due to the shape of the mass function close to the 
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flux limit of the different surveys and redshifts considered 
here. 



6 DISCUSSION 

6.1 Strong self-similarity 

The examination of the slopes of the Lx — kT relation in 
different sub-populations of this large cluster sample indi- 
cated that the most relaxed systems, or those hosting the 
strongest cool cores (definitions with significant overlap), ex- 
hibit a low scatter Lx — kT relation with a slope of ~ 2, in 
agreement with the self-similar model, once their core re- 
gions have been excised. Such a self-similar correlation has 
generally not been found in previous studies of the Lx — kT 
relation, which is likely due our large sample size, and to 
the fact that most previous studies have not distinguished 
CC and NCC s ubsamples once a core correction has been 
made. However, IPratt et al.l l|2009l . hereafter P09) did make 
such a distinction in their recent XMM-Newton study of the 
Lx — kT relation of the REXCESS sample, and used a simi- 
lar methodology to us allowing straightforward comparisons 
to be made. 

In the P09 study, after core exclusion, the Lx — kT 
slopes of the relaxed/CC subsamples were shallower than 
the unrelaxed/NCC subsamples, but (in contrast to our re- 
sults) were still significantly steeper than the self-similar 
slope of 2. The methods used to define the CC/NCC and 
relaxed/unrelaxed classes in P09 differ slightly from those 
employed here. P09 separated CC and NCC clusters using 
a measurement of the central density, and classed ~ 30% of 
their clusters as CC (a similar fraction to us). To classify 
relaxed and unrelaxed clusters, P09 used the same centroid 
shift measurement as us (albeit with a slightly different im- 
plementation), but used a value of (w) = O.OI-R500 to seg- 
regate the relaxed/unrelaxed subsets. This less strict defini- 
tion of a relaxed cluster than our (w) > O.OO6-R500, resulted 
in a significantly larger fraction of the P09 sample (~ 60%) 
being classed as relaxed. 

Subsample definition and instrumentation notwith- 
standing, the key difference between the P09 and our sam- 
ples is the temperature range of systems covered. Half of the 
P09 CC clusters are cooler than 3.5 keV (our coolest CC 
system) . This is made clear in Figure 1121 which plots the 
relaxed and CC P09 clusters on the corresponding Lx — kT 
relations of our sample. At kT> 3.5 keV the P09 data agree 
quite well with our Lx — kT relations, with their essentially 
self-similar slopes (particularly for the CC selection). This 
suggests that the steeper P09 Lx — kT slope is driven by 
their inclusion of lower mass systems, and that the self- 
similarity of relaxed/CC clusters found here is limited to 
mor e massive systems (k T > 3.5 keV). 

ICroston et al.l (|2008T ) presented the gas density profiles 
of the REXCESS clusters, allowing us to investigate if this 
possible similarity breaking of the Lx — kT relation below 
3.5 keV is manifested in the gas density profiles. In Fig. \T3\ 
we show the scaled gas density profiles of the CC REX- 
CESS clusters, separated into hot (kT > 3.5 keV) and cool 
(kT < 3.5 keV) subsets. Also indicated is the mean profile 
of the relaxed clusters from our sample (as plotted in Fig[6} . 
The hot REXCESS profiles agree fairly well with our mean 
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Figure 12. Top: the core-excised Lx — kT relation for the CC 
clusters in the current sample is plotted along with the CC clus- 
ters of P09. Bottom: the core-excised Lx — kT relation for the 
relaxed clusters in the current sample is plotted along with the 
relaxed clusters of P09. 

relaxed profile (which are also all kT > 3.5 keV clusters), 
but the cool REXCESS CC clusters have flatter profiles, 
with lower gas densities than the other cluster profiles out 
to ~ 0.7i?5oo- This suggests that the self-similarity observed 
in the gas density profiles of the relaxed clusters in our sam- 
ple is valid only for hotter (kT > 3.5 keV) systems. This 
is consistent with the proposed steepening of the Lx — kT 
relation below 3.5 keV. 

The scaled gas density profiles of our large sample also 
provide information on where self-similarity breaks down. 
For the unrelaxed/NCC clusters, these profiles showed that 
the similarity index was > out to a cluster-centric radius 
of « 0.7/?5oo- Beyond this radius, the gas density profiles of 
the unrelaxed/NCC clusters appear self-similar. These re- 
su lts broadly ag r ee wi th the XMM-Newton measurements 
of lCroston et al. (|2008l . although in that work, clusters were 
not separated into relaxed/unrelaxed subsamples). In con- 
trast, the similarity index of the relaxed/CC clusters is con- 
sistent with zero at all radii probed (Fig. [9|. This would 
seem to imply that excluding the cores should not be re- 
quired to recover a self-similar Lx — kT relation. This may 
be the case, but as illustrated in Figures [8] while there is 
no strong temperature dependence of the scaled gas den- 
sity in the cores, there is departure from self-similarity in 
the form of large cluster-to-cluster variation in scaled gas 
density. This is indicative of different cool core strengths 
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Figure 13. Scal ed gas density profil es of the REXCESS CC clus- 
ters (taken from lCroston et al.ll2008h . The REXCESS clusters are 
divided into hot and cool subsets as indicated. The hatched re- 
gion indicated the standard deviation of the scaled gas density 
profiles of the relaxed clusters in our sample. 



and is likely responsible for the non self-similarity of the 
core-included Lx - kT relation. Finally, the REXCESS CC 
profiles in Fig. 1 131 show that in cool (kT < 3.5 keV) CC sys- 
tems, the core gas density is lower than in hotter systems, 
suggesting that a non-zero similarity index would be mea- 
sured at radii smaller than ~ 0.7-Rsoo for the relaxed/CC 
clusters in our sample if it included cooler systems. 

Outside a radius of ~ 0.15i?soo, the similarity index of 
relaxed/CC clusters is close to zero and the cluster-to-cluster 
variation in scaled gas density is low, indicating that the 
ICM in these regions is self-similar. This is consistent with 
the self-similar slope of the core-excised Lx — kT relation of 
the relaxed/CC clusters. 

Our results can a l so be usefully compared with the 
study by iMantz et al.l (|2010ah who found that when core 
regio ns were excis e d, the Lx — kT relation slope was 2.70 ± 
0.20. IMantz et ail (|2010sJ) did not separate relaxed and dis- 
turbed or CC/NCC clusters in their analysis, but their slope 
is in excellent agreement with the slope we fin d for all clus- 
ters wit h cores excluded (2.72 ±0.18). Usefully. IMantz et alj 
l|2010al ) were able to include full treatment of selection biases 
in their analysis, and showed that for core-excised luminosi- 
ties, selection biases did not significantly impact the mea- 
sured scaling relations, so the self-similar slope we find for 
the core-excised Lx — kT relation of the relaxed/CC clusters 
in our sample is not likely to be affected by selection bias. 

We thus conclude that the ICM of massive (kT> 
3.5 keV), relaxed (or CC) clusters obeys strong self- 
similarity outside the core (r < 0.15i?soo) regions. 



6.2 Implications for similarity-breaking models 

iPratt et al.l l|201Ch argued that variations of the gas con- 
tent of clusters with mass and radius are at the root of the 
observed departures from self-similarity of cluster entropy 
profiles. They proposed a mechanism whereby gentle heat- 
ing from AGN, coupled with merger related mixing is re- 
sponsible for a redistribution of gas, leading to the observed 
entropy profiles, and also the suppression of luminosity in 
low mass systems, giving rise to steepening in the Lx — kT 



relation. Our results are consistent with this picture, but add 
the requirement that the most massive relaxed/CC systems 
are essentially self similar outside t he central . 15fl5q o. The 
entropy of the CC clusters in the IPratt et al.l (|201Ch sam- 
ple shows evidence for excess above the self-similar expecta- 
tions at radii as large as ~ 0.7i?soo, however, this is driven 
by the lower mass systems. Clusters with kT> 3.5 keV are 
in reasonable agreement with the self-similar entropy model 
be yond ~ 0.15/?5oq (roughly O.I-R200, as plotted in figure 3 



Pratt et alj|20ich . 



Our results on the details of the Lx — kT relation and 
density structure of the ICM present new information with 
which to refine models of ICM heating and feedback. In par- 
ticular, by identifying a fully self-similar regime in the clus- 
ter parameter space, we have a baseline against which to 
measure the impact of feedback. 

Let us first consider the case of relaxed/CC clus- 
ters. With the core regions excluded, we find a self-similar 
Lx — kT relation, which appears (when combined with the 
P09 data) to steepen below ~ 3.5 keV. The density pro- 
files of the kT > 3.5 keV relaxed/CC systems are also self- 
similar beyond the core, with low dispersion and similarity 
index close to zero, but again the cooler REXCESS systems 
show similarity breaking, with suppressed densities relative 
to the hotter systems. Our results thus suggest that for re- 
laxed systems above ~ 3.5 keV, the effects of feedback are 
negligible outside the core of 0.15i?5oo- 

A possible interpretation of this is that feedback occurs 
in the cores of all of the relaxed/CC systems, providing some 
fixed level of energy input beyond the balancing of the gas 
cooling. The heated gas will have raised entropy, and will 
rise, expand and cool, conserving entropy, until it reaches 
a radius of gas with the same entropy level. For more mas- 
sive syste ms, this will occu r at a smaller radius (see e.g. 
figure 1 in IPratt et alj|2010h . Thus for a mass-independent 
amount of heating, there will be a mass threshold above 
which the effects of feedback are confined to the core re- 
gions. For a 3.5 keV system, the en tropy level at 0.157?5oo is 
~ 200 keV cm 3 (|Pratt et alJ feoiO). so an entropy increase 
somewhat lower than this would not redistribute gas much 
beyond the core. In the absence of mergers, the new entropy 
structure will remain intact until "eroded" by radiative cool- 
ing, but periodic core feedback would maintain a quasi-static 
ICM configuration. 

For unrelaxed/NCC clusters, the population does not 
appear to be self-similar at any mass, even with the cores 
excluded. This is manifested in the steep Lx — kT relation 
which crosses that of the relaxed/CC clusters at ~ 6 keV, 
and also in the gas density profiles, which show tempera- 
ture dependence out to large scaled radii. The steepening 
of the Lx — kT relation below 6 keV can be explained in 
the same terms as for the < 3.5 keV relaxed/CC clusters 
if we posit that core feedback occurs in all clusters (not 
just relaxed/CC) clusters. This is consis tent with quasar 
mode heat input occurring in all clusters (IShort fc Thomas! 
2009). We then require that the effects of feedback be no- 
ticeable outside the core for more massive systems (up to 
kT ~ 6 keV). This requires a larger increase in entropy, 
as the ambient entropy level at O.I5-R500 in a 6 keV sys- 
tem is ~ 300 keV cm 3 l|Pratt et all 120101 ). This is achiev- 
able in principal without additional feedback energy, via 
the amplification by merger shocks of entropy input in clus- 
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ter cores. This amplification is particularly effective if the 
original entropy increase was driven by lowering of the gas 
density rath er than increasing its temperature (this e ffect is 
discussed in lPonman. Sanderson fc Finoguenovll2003l . in the 
context of the accretion of material during cluster forma- 
tion). With the entropy further increased by merger shocks, 
the gas may then be able to settle to radii well beyond 
the cluster cor e . Thi s model is similar to that proposed 
by iPratt et al.l (|201Cf ). but with an emphasis on the role 
of mergers to amplify the entropy increase due to central 
energy input. 

A challenge remains if we wish to explain all aspects 
of the observed departures from strong self-similarity in our 
sample with this central heating and merger shock mecha- 
nism: the loci of the > 6 keV unrelaxed/NCC clusters lie 
above the self-similar Lx — kT relation. We thus speculate 
that while the combination of central heating and merger 
shocks may act to raise gas entropy beyond 0.15-Rsoo in 
kT < 6 keV systems, at higher masses the dominant ef- 
fect of mergers is to enhance Lx and/or suppress kT. We 
note that numerical simulations o f clusters do not unam- 
biguou sly support this hypothesis. iRowlev. Thomas fe Kavl 
Hooil) found that merging clusters tend to move along the 
Lx — kT relation, though it was observed that Lx increases 
and kT decreases when the cores of two systems merge. A 
similar effect was also present in the merger simulations of 
iPoole et al.1 l|2007h , in which massive > 6 keV clusters were 
shown to have enhanced Lx relative to kT when the cores 
are at their closest initial approach, before settling down 
onto the Lx — kT relation over a period of w 4 — 5 Gyr. 
This could explain the observed position of the unrelaxed 
clusters above the self-similar relation at the high-kT end, 
particularly combined with a selection bias to preferentially 
detect and observe clusters that are at the more luminous 
stages of their mergers (although we note that at a close core 
passage, such clusters may not be identified as unrelaxed). 

By cont r ast, in the cosmological simulations of 
lHartlev et al.l |2008), mergers were found to push clusters 
along, but slightly below the mean Lx — kT relation, lead- 
ing to a curve towards a flatter Lx — kT relation for dis- 
turbed clusters at the high kT end. This appears to contrast 
with our results, but the comparison is not straightforward, 
as the simulations are necessarily dominated by lower mass 
(kT< 3 keV) objects due to the limited volume simulated, 
and the classification of disturbed clusters is based on time 
since the last major merger, rather than a morphological 
measurement. For the same reasons, the results of these dif- 
ferent computational studies are not necessarily at odds. 

6.3 Heating or Lx suppression? 

Given the loci of lower mass systems below the self-similar 
Lx — kT relation, it is interesting to consider whether those 
clusters occupy that position due to their temperatures be- 
ing enhanced relative to their virial temperature, or their lu- 
minosities being suppressed due to restructuring of the ICM. 
In some sense, this is just a question of timescales, as heat- 
ing of a parcel of gas will eventually (on a timescale of order 
the s ound-crossing time: a few 10 8 years; iMathews fc Guol 
l201ll ) have risen, expanding and cooling, to its new adi- 
abat, resulting in a restructuring of the ICM. Most theo- 
retical work focuses on this longer term result of the re- 



structured ICM, but some observational studies have ar- 
gued that for clusters hosting radio-loud AGN, the steepen- 
ing of the Lx — k T relation is due to temperature enhance- 
ment of the ICM llCroston. Hardcastle fc Birkinshawll2005l : 
iMagliocchetti fc Briiggenl 120071 ). 

In the absence of a coherent set of radio data for these 
clusters, we can appeal to their gas density profiles to inves- 
tigate whether the steepening of the Lx — kT relation is due 
to enhanced kT or suppressed Lx of the ICM. Fig. H3lshows 
that the ICM structure of the lower mass CC clusters is 
different from the self-similar higher mass CC profiles. The 
ICM density in the lower mass CC clusters is lower than or 
similar to the self-similar profiles, suggesting that the lumi- 
nosity is suppressed in these systems by the removal of gas 
from the central regions out to R500 and beyond. 

For unrelaxed clusters, Fig. |S] shows that their ICM 
structure is flatter than that of the self-similar relaxed 
clusters. Indeed, the gas density is higher for unrelaxed 
clusters at scaled radii J>0.4i?5oo- A similar effect is 
seen in the scaled gas density profiles of the Planck- 
selected sample of largely unrelaxed clusters, when com- 
par ed with the mean profile of the full REXCESS sam- 
ple l|Planck Collaboration et al.ll201lT ). Indeed, the same ex- 
cess of gas density at large radii in NCC clusters is also 
present (albeit more weakly) in the REXCESS profiles of 
|Croston et al.ll2008h if they are separated into CC and NCC 
subsets. Given that the gas density is higher than the self- 
similar profiles at large scaled radii, it is not immediately 
clear if the net effect would be a suppression or increase of 
Lx- This can be crudely estimated by defining a "psuedo- 
Lx" as the integral of the square of the scaled mean den- 
sity profile over the spherical shell from [0.15 — lJ-Rgoo- This 
calculation shows that the unrelaxed clusters have a mean 
luminosity ~ 15% lower than the relaxed clusters, so that 
the net effect of the observed ICM structural differences is 
suppression of Lx. 

This picture changes slightly if the unrelaxed clusters 
are split into hot (kT > 6 keV) and cool (kT < 6 keV) sub- 
sets. The mean profiles of these subsets are shown in Fig. 
1141 The temperature dependence of the unrelaxed profiles is 
apparent, with cooler systems having a flatter mean profile, 
which crosses the mean relaxed profile at a larger scaled ra- 
dius than the hotter unrelaxed systems. This is consistent 
with the picture of core feedback resulting in the removal of 
gas to larger radii in lower mass systems. The effect on the 
system luminosity is also different. For the hotter unrelaxed 
clusters, the psuedo-Lx is ~ 5% higher than the relaxed 
clusters, while for the cooler unrelaxed clusters it is ~ 20% 
lower. This is qualitatively consistent with the crossover of 
the relaxed/CC and unrelaxed/NCC Lx — kT relations at 
6 keV. For a quantitative comparison, the ratio of the lumi- 
nosity predicted by the RCC and NRCC relations (Fig. [5} 
was computed for the temperature of each unrelaxed clus- 
ter. For the kT < 6 keV clusters, the mean ratio was 0.74, 
indicating that the 20% reduction in Lx suggested by the 
gas density profiles is sufficient to explain the steeper unre- 
laxed/NCC relation, without enhancement of the ICM tem- 
perature. For the kT > 6 keV unrelaxed clusters, the mean 
ratio was 1.30, significantly larger than the 5% increase in 
Lx predicted by the mean density profiles. 
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Figure 14. Average scaled gas density profiles of the relaxed 
and unrelaxed subsamples, with unrelaxed clusters subdivided 
into hot (kT > 6 keV) and cool (kT < 6 keV) subsets. 
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Figure 15. Histogram of F CO re values for clusters at z < 0.5 
(lines) and z > 0.5 (filled bars). Each histogram is normalised to 
unity. 



6.4 Weak self-similarity 

Our investigation of the evolution of the Lx — kT relation 
suggests at first glance, a set of clusters at intermediate 
redshifts (0.2 < z < 0.6) with luminosities significantly 
higher than predicted by self-similar evolution of the low- 
re dshift relation (Figure I10[) . This is similar to the results 
of lBranchesi et al. (2007), who argued in favour of non self- 
similar evolution. Our simple simulations of a self-similar 
population of clusters subjected to selection functions that 
plausibly approximate those of our heterogeneously selected 
sample indicated that the deviations from self-similarity 
can reasonably be explained by selection biases. As demon- 
strated by the simulations, the combination of intrinsic scat- 
ter <Jl\m, and the slope of the mass function at the mass cor- 
responding to the survey flux limit give differing amounts of 
bias for different flux limit and redshift combinations. We 
thus find no evidence to reject the self-similar description 
of evolution in the Lx — kT relation. This conclusion is in 
agreement with other analyses based on statistically com- 
plete samples where full tr eatment of selection effects is pos- 
sible dPacaud et all (|2007l . 29 clusters at z <; 1) iMantz et al] 
(|2010al . 238 clusters at z<,0.5)) which found no strong 
evidence for departures from self-similar evolution in the 
Lx — kT or Lx — M relation (althoug h we note that in a 
similar analysis, IVikhlinin et al. 2009a, found an indication 
of evolution slightly weaker than the self-similar prediction 
for the Lx — M relation). 

Our simple simulations assumed ctl\m was constant 
with r edshift, but in fact this may not be the case. iMaughanl 
l|2007l ) showed that <tl\m was significantly lower for the clus- 
ters in the current sample at z > 0.5 than those at z < 0.5. 
This is consistent with evidence that the fraction of clus- 
ters with cool cores (the domina nt contributor to Ql\m ) 
decreases strongly above z w 0.5 (Vik hlinin et al.l l2006al ). 
Evolution in (Tl\m could introduce further bias in measure- 
ments of the Lx — kT relation evolution. This evolution of 
the CC clusters was investigated in our sample, and Figure 
1151 shows normalised histograms of F core for z > 0.5 and 
2 < 0.5 clusters. 

Fig. ll5l shows a clear absence of strong CC (F COIC > 0.5) 
clusters at z > 0.5 in the cu rrent samp l e. A similar trend 
is seen in the simulations of iKav et al.l |2007l ). who found 



a corresponding red uction in ctl,\t with red shift. Our re- 
sults agree with the IVikhlinin et al. (2006a) findings, and 
we also found similar histograms when using the cuspiness 
CC proxy. We also showed previously that the fraction of 
relaxed clust ers decreased sig nificantly at z > 0.5 in the cur- 
rent sample (Maughan 2007). However, while these results 
provide strong evidence for evolution in the CC population 
in our sample, we remind the reader once more of the hetero- 
geneous nature of our sample and urge caution in extrapolat- 
ing these results to the general cluster population (see e.g . 
ISanderson. O'Sullivan fc Ponmanl 120091 ; iMittal et ail l201ll . 
for discussions of cool core properties in representative sam- 
ples). Indeed, even in statistically complete samples, the 
fraction of CC clusters at high redshift is subject of some 
debate. ISantos et al.l (|2010l ) recently showed that the cool 
core fractions in three high-redshift ROSAT-derived surveys 
differed significantly, most likely due to their different clus- 
ter detection algorithms (high-z cool core clusters could be 
classed as point sources and excluded from samples) . In two 
of the surveys considered, the distribution of CC strengths 
(measured with a method similar to our Feme proxy) was 
similar at low and high [z > 0.6) redshift. This would imply 
weaker (or zero) evolution in (?i,\m than found in our non- 
representative sample. The uncertainty on the evolution in 
&l\ M further complicates attempts to fully include selection 
bias corrections in studies of the evolution i n the Lx — kT 
relation, although note that I Vikhlinin et all l|2009al ) showed 
that estimates of a jam from flux-limited samples are unbi- 
ased, making bias-corrected studies feasible. 

The redshift evolution of the relaxed/CC fraction in 
our sample means that the relaxed/CC and unrelaxed/NCC 
subsets cover different redshift ranges with the relaxed/CC 
subsets essentially limited to z < 0.5. In order to test 
whether this difference contributed to the shallower Lx — kT 
slopes found for the relaxed/CC clusters, the Lx — kT rela- 
tions of the unrelaxed/NCC clusters were fit for z < 0.5 and 
z > 0.5 clusters separately. There slopes of both the high 
and low redshift unrelaxed/NCC clusters agreed very well 
with the fit to the full unrelaxed/NCC subsample, and there 
was no evidence for a change of slope with redshift. 
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7 CONCLUSIONS 

We have used a sample of 114 clusters observed with Chan- 
dra ACIS-I across a wide baseline in temperature (2 < kT < 
16 keV) and redshift (0.1 < z < 1.3) to study the self simi- 
larity of the cluster population in terms of its mass scaling 
and redshift evolution. Our main conclusions on the mass 
scaling of the clusters are as follows: 

• The ICM of massive (kT > 3.5 keV), relaxed/CC 
galaxy clusters is self-similar outside of the central 0.15-Rsoo- 
In this regime, strong self-similarity is obeyed, manifested in 
an Lx — kT relation with a slope of » 2, and ICM density 
profiles with low dispersion and no temperature dependence. 

• By comparing our data with measurements of the REX- 
CESS sample, which extends to lower masses, we find that 
the self similarity of the relaxed/CC clusters breaks below 
~ 3.5 keV, manifested by a steepening Lx — kT relation and 
flatter density profiles. This implies that the impact of cen- 
tral heating extends beyond the core in these lower mass 
systems. 

• Unrelaxed/NCC clusters are not self-similar; their 
Lx — kT relation has a steeper than self-similar slope, and 
their ICM density profiles are temperature-dependent out 
to » 0.77?500- 

• The steeper unrelaxed/NCC core-excised Lx — kT re- 
lation crosses the self-similar relaxed/CC relation at around 
6 keV. Below this temperature unrelaxed/NCC clusters ap- 
pear to be less luminous than their relaxed/CC counter- 
parts, while above this temperature they appear to be cooler 
and/or more luminous. 

• These results are consistent with similarity breaking in 
clusters being due to central feedback, the effects of which 
extend beyond the central 0.15J?soo in low mass ( <, 3.5 keV) 
systems. We suggest that merger shocks act to amplify the 
entropy increase from this feedback, allowing its effect to be 
felt beyond the core in more massive ( <, 6 keV) unrelaxed 
systems. In the most massive systems, any heating is limited 
to the core regions and the dominant effect of mergers is to 
raise the luminosity of the ICM or lower its temperature 
relative to the self-similar Lx — kT relation. 

Our investigation of the evolution of the Lx — kT rela- 
tion gave the following main results: 

• The evolution of the Lx — kT relation in our sample is 
inconsistent at face value with the self-similar model, but 
these differences are plausibly explained by a reasonable 
model of the varying selection biases due to the different 
selection functions in place across our heterogeneous sam- 
ple. 

• Our bias modelling assumed the intrinsic scatter in the 
Lx — M relation (o L \ M ) is redshift independent, however, in 
our sample the fraction of cool core clusters is much lower 
at 2 > 0.5, consistent with th e smaller < jl\m w e found at 
z > 0.5 for the same clusters in lMaughanl (|2007T ). Such vari- 
ation in o"i,|m could further bias evolution measurements, 
but should be applied to the general population with cau- 
tion due to the non-representative nature of our sample. 

• We would thus argue that for the core-excised Lx — kT 
relation, self-similar evolution is obeyed to at least a first or- 
der approximation. This suggests that the balance of heat- 
ing, cooling and mergers has remained roughly constant 
since z ~ 1. Stronger constraints on non self-similar be- 



haviour require statistically complete samples at all red- 
shifts, with a full treatment of the selection function, in- 
cluding the sensitivity of the cluster detection algorithm to 
cool core clusters at high redshift. 

C urrent X-ray surveys such as XXL dPierre et al.l201Ch 
XCS i|Llovd-Davies et all l2O10h . and XDCP (|Fassbenderl 
2008) are beginning to produce suitable samples to 
study cluster evolution on a stronger sta ti stical foun- 
dation, extending the work of iMantz et al.l (|2010aT ) and 
IVikhlinin et all |2009af ) to z > 1. Further insight will be 
gained by studies of samples selected independently of the 
X-ray emission such as optical, weak lensing, and to some 
extent the Sunyaev-Zel'dovich effect (SZE; although the X- 
ray properties of SZE selected clusters may retain some bias 
due to the dependence of both on the physical properties of 
the ICM). 
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